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ABSTRACT 

Aims. By means of hydrodynamical models we do the first investigations of how the properties of planetary nebulae are affected by 
their metal content and what can be learned from spatially unresolved spectrograms of planetary nebulae in distant stellar systems. 
Methods. We computed a new series of ID radiation-hydrodynamics planetary nebulae model sequences with central stars of 
0.595 Mq surrounded by initial envelope structures that differ only by their metal content. At selected phases along the evolutionary 
path, the hydrodynamic terms were switched off, allowing the models to relax for fixed radial structure and radiation field into their 
equilibrium state with respect to energy and ionisation. The analyses of the line spectra emitted from both the dynamical and static 
models enabled us to systematically study the influence of hydrodynamics as a function of metallicity and evolution. We also recom- 
puted selected sequences already used in previous publications, but now with different metal abundances. These sequences were used 
to study the expansion properties of planetary nebulae close to the bright cut-off of the planetary nebula luminosity function. 
Results. Our simulations show that the metal content strongly influences the expansion of planetary nebulae: the lower the metal con- 
tent, the weaker the pressure of the stellar wind bubble, but the faster the expansion of the outer shell because of the higher electron 
temperature. This is in variance with the predictions of the interacting-stellar-winds model (or its variants) according to which only 
the central-star wind is thought to be responsible for driving the expansion of a planetary nebula. Metal-poor objects around slowly 
evolving central stars become very dilute and are prone to depart from thermal equilibrium because then adiabatic expansion con- 
tributes to gas cooling. We find indications that photoheating and line cooling are not fully balanced in the evolved planetary nebulae 
of the Galactic halo. Expansion rates based on widths of volume-integrated line profiles computed from our radiation-hydrodynamics 
models compare very well with observations of distant stellar system. Objects close to the bright cut-off of the planetary nebula 
luminosity function consist of rather massive central stars (> 0.6 Mq) with optically thick (or nearly thick) nebular shells. The half- 
width-half-maximum velocity during this bright phase is virtually independent of metallicity, as observed, but somewhat depends on 
the final AGB-wind parameters. 

Conclusions. The observed expansion properties of planetary nebulae in distant stellar systems with different metallicities are ex- 
plained very well by our ID radiation-hydrodynamics models. This result demonstrates convincingly that the formation and accel- 
eration of a planetary nebula occurs mainly because of ionisation and heating of the circumstellar matter by the stellar radiation 
field, and that the pressure exerted by the shocked stellar wind is less important. Determinations of nebular abundances by means of 
photoionisation modelling may become problematic for those cases where expansion cooling must be considered. 

Key words, hydrodynamics - planetary nebulae: general - planetary nebulae: individual (M2-29, NGC 1 360, NGC 4361 , NGC 7027, 
PNG135.9+55.9) - stars: AGB and post-AGB) 



1. Introduction 

The use of planetary nebulae (PNe) has become an important 
tool for investigating the properties of unresolved stellar popu- 
lations in galaxies. Applications are kinematic studies for prob- 
ing the gravitational potential, luminosity functions for distance 
estimates, and, most importantly, abundance determinations by 
means of plasma diagnostic or photoionisation modelling. 

It has not been thoroughly investigated to date whether re- 
sults based on PNe are generally trustworthy. For instance, by 
applying standard plasma diagnostics, it is implicitly assumed 
that a PN is a homogeneous object in which all physical pro- 
cesses are in equilibrium. In reality, however, a PN is a highly 
structured expanding gas shell, subject to the dynamical pres- 
sure of a fast central-star wind and the thermal pressure of the 
photo-heated gas. Only if radiation processes dominate heating 
and cooling of the gas hydrodynamics can be neglected. The 
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questions whether and by how much PNe properties, and espe- 
cially thermal equilibrium, depend on their metal content, has 
not been posed yet. 

The recent progress in modelling the evolution of PNe by 
means of radiation-hydrodynamics opens the possibility of a bet- 
ter appreciation of time-dependent effects. Marten ( 1993a, 1995) 
demonstrated that non-equilibrium may become important for 
rapidly expanding, low-density nebular shells, and also for cases 
with quickly evolving central stars. Well-known examples are 
the hot PN haloes that are extremely out of thermal and ionisa- 
tion equilibrium ( |Marten| 1993b] |Sandin et al.|2008 1. 

Later |Perinotto et aT] ( [T99 8 ) used hydrodynamical models to 
check the influence of hydrodynamics on the abundance deter- 
mination by plasma diagnostics. The authors found no signifi- 
cant non-equilibrium effects, but note that large abundance er- 
rors may occur in more advanced evolutionary stages because of 
obviously inappropriate ionisation correction factors. However, 
since the number of model sequences used by the authors was 
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only small and restricted to a metallicity typical for Galactic 
disk objects, the case concerning the reliability of abundance 
determinations by plasma diagnostics is far from being set- 
tled. In a forthcoming paper we will discuss the reliability of 
ionisation correction factors presently used in more detail by 
means of the hydrodynamical simulations presented in this paper 
(Schonberner et al., in prep.). 

While density and velocity structures of planetary nebu- 
lae are clearly the result of the dynamics, the case is differ- 
ent for the electron temperature. Deviations from the equilib- 
rium value may occur by dynamical effects, i.e. by local gas 
compression/expansion, and/or by radiative effects, i.e. by heat- 
ing/cooling and ionisation/recombination. The radiative effects 
occur either by rapid stellar evolution and/or by changes of 
the optical depth (Marten 1995 1. Since temperature deviations 



from its equilibrium value cannot unequivocally be attributed to 
only one of these physical processes, we simply speak of "non- 
equilibrium" effects. 

To date, existing radiation-hydrodynamics simulations used 
the typical chemical composition of Galactic disk PNe. In the 
context of the grown interest in extragalactic studies, also of PNe 
populations, it is of great importance to investigate how the evo- 
lutionary properties of PNe will change if they originate from a 
metal-poor population. 

The metal content is expected to influence the PN evolution 
twofold: Firstly, the cooling efficiency of the gas decreases with 
metallicity, leading to a higher electron temperature. Secondly, 
the strength of the central-star wind decreases with metallicity 
because of the reduced line opacity which drives the wind. PNe 
in stellar populations with low metal content are thus expected 
to expand faster than their Galactic disk counterparts since the 
expansion speed is roughly proportional to the sound speed, i.e. 



to the square root of the electron temperature (Schonberner et 
[aL] [200 5a hereafter Paper II). Together with a reduced wind 
power, metal poor PNe will be substantially more dilute than 
their Galactic counterparts and may thus quickly run into non- 
equilibrium since the line cooling efficiency of the gas depends 
on the density squared and may become less important for the 
nebular energy budget than expansion cooling. The radiation 
field of the central star, and to some extend also its luminosity, 
are in principal metal dependent, too. Both effects are usually 
not considered and are neglected here for simplicity. 

Already from the sample of planetary nebulae of the Galactic 
halo one can get observational hints that low-metallicity PNe are 
prone to deviate from thermal equilibrium. This sample com- 
prises objects in different evolutionary stages with metallicities 
below solar at various degrees. The study of halo PNe appears 
very promising since they can be used as proxies for extragalac- 
tic nebulae which are generally too faint and too small for spa- 
tially resolved investigations with high spectroscopic resolution. 

|Howard et al.| ( |T997] > analysed the then known 9 PNe of the 
Galactic halo in a very homogeneous way by collecting all the 
available line data from the literature and by fitting the spectra 
with photoionisation models such that "... a reasonable match 
between model and observations is achieved". Closer inspection 
of the results listed in |Howard et al.| ( [l997[ Table 2 therein) re- 
veals clearly that the objects with the highest electron temper- 
atures show a remarkable inconsistency between the observed 
and predicted ratios of the temperature sensitive [O m] lines, i.e. 
for these cases the procedure of a general line-strength fitting 
from the ultraviolet to the optical spectral range did not lead to a 
consistent determination of a (mean) electron temperature! 

The case is illustrated in Fig.[T] where we plotted the ra- 
tio of observed and modelled [O m] lines against observed val- 
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Fig. 1. Ratio of observed and modelled strengths of relevant 
[Om] lines (indicated in the legend) vs. the temperature sensi- 
tive quantity Rom - (/5007 + ^4959)/A363 for the 9 Galactic Halo 
PNe investigated by Howard et al. ( 1997] >. 



ues of the temperature indicator Ro m = (/5007 + ^4959)//4363- We 
see that for lower electron temperatures (larger Rom) the pho- 
toionisation models are able to reproduce the observed [Om] 
line ratios rather well within the 10 % error level. However, for 
Rom 55 60 (or T e ^ 16 000 K) discrepancies between observa- 
tions and predictions become very large. 

These differences are especially seen for the temperature 
sensitive A 4363 A line which is off by up to 30-40% for 
Rom ^ 60^]We interprete this behaviour at low Rom as a clear 
indication for the failure of photoionisation modelling because 
of a possible breakdown of thermal equilibrium at low metallic- 
ities. 

Concerning the expansion properties of PNe and their ex- 
pected dependence on metallicity, the present situation is also 
unclear. Richer (2006) studied the expansion behaviour of PNe 
in Local Group galaxies and could not find any significant de- 
pendence on the metal content: the brightest PNe always have, 
on the average, a HWH M velocity (in the [O in] /I5007 A line) 
of ^18 kms~' (see also |Richer||2007j ). A similarly low value 
(16.5 kms 1 ) has recently been found for the brightest PNe in 
the Virgo cluster by Arnaboldi et al. (2008). This finding is puz- 
zling since theory predicts, as outlined above, a clear correlation 
of the PN's speed of expansion with metallicity. 

Because of the importance of planetary nebulae for extra- 
galactic studies, we decided to investigate how planetary nebulae 
evolution depends on metallicity. In particular, we are interested 
in the expansion velocities of the brightest PNe and how these 
correlate with the metallicity of the parent stellar population, and 
in the question under which conditions standard photoionisation 
models may fail. 

We start with describing the setup of our ID radiation- 
hydrodynamic models in the next section, and continue in Sect. [3] 



1 We have no explanation why the modelled A 4363 A line is some- 
times stronger and sometimes weaker than the observed line. It is prob- 
ably the outcome of the fitting procedure. 
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Table 1. Distribution of the chemical abundances typical for 
Galactic disk PNe. The abundances, e,, are given in the usual 
manner as (logarithmic) number fractions relative to hydrogen, 
i.e. e ; = log(w,/n H ) + 12. 
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with a detailed presentation of our results. In Sect. |4]we use our 
models to interprete the expansion velocities observed in distant 
stellar systems in terms of metallicity and evolutionary stage. 
The paper is concluded with Sects. |5]and|6] Part of the results 
presented in Sect.[3]can be found in Schonberner et al. ( 2005c| l. 



2. The models 

Our method of modelling the evolution of planetary nebulae 
by means of ID radiation-hydrodynamics has been documented 
in previous publications and shall not be repeated here (see 
Peri notto et al.||2004] and references therein; henceforth called 
Paper I). A very detailed comparison of our models with ob- 
served nebular structures can be found in Steffen & Schonberner] 
(2006). We add only that all hydrodynamical sequences used in 
this work are based on a new parallelised version of the code 
which also includes thermal conduction as described in detail 
by Steffen et al. (2008). Already existing sequences have been 
recomputed with this new version in order to achieve a homo- 
geneous set of evolutionary sequences. We emphasise that, once 
the initial envelope structure and the central star, whose radiation 
field and wind represent the time-dependent inner boundary con- 
ditions, have been selected, the whole evolution is consistently 
determined by the hydrodynamics. We note that all the central- 
star models used in this work burn hydrogen, which means that 
the case of nebulae around hydrogen-deficient, i.e. Wolf-Rayet, 
central stars is thus not considered here. 

Concerning the elements to be considered in the nebular en- 
velopes, we included, next to hydrogen and helium, only the 
most important coolants, i.e. carbon, nitrogen, oxygen, neon, 
sulphur, chlorine, and argon. For each individual element up to 
12 ionisation stages are taken into account. We emphasise here 
that ionisation, recombination, heating and cooling are treated 
fully time dependently, and the (radiative) cooling function for 
each volume element is composed of the contributions of all in- 
dividual ions (see also Marte n~& Szczerbap 997). Our reference 
abundance distribution for the Galactic disk PNe, Zgd, is listed 
in Table Q] These abundances are the same as those used in our 
previous hydrodynamical computations and are, apart from car- 
bon and nitrogen, very close to the most recent solar values (see, 
e.g., |Asplundet al.|2009| >. 

We considered in general six cases with scaled Galactic disk 
abundances distributions, Zqd, of C, N, O, Ne, S, CI, and Ar, 
covering the range from Z = 3Zgd to Z = Zgd/100 in steps 
of 0.5 dex. For simplicity, we did not consider metallicity- 
dependent variations of abundance ratios in this pilot study. 
The range of metal contents considered here covers the ob- 
served degrees of metallicities in galaxies: While giant galax- 
ies have metal contents within a factor of 3 of the solar one, 
dwarf galaxies have less metals, usually about 1/10 of the so- 
lar case. The metal-poorer sequences are thought to demonstrate 
non-equilibrium effects more clearly. They are useful for inter- 
preting objects with very low metallicity, like, for instance, the 
Galactic halo object PN G135.9+55.9 which has a mean metal 




Fig. 2. Radial run of ion density (thick) and gas velocity (thin) 
for an initial nebular configuration with a = 3 density profile at 
age zero. The outer boundary is at r — 2.8 x 10 18 cm (outside the 
graph). The stellar parameters are given at the panel's top and 
refer to the 0.595 Mq post-AGB track shown in Fig. [3] Density 
and velocity at r — 1.0 x 10 16 cm correspond to an AGB mass- 
loss rate of 1.3 x 10~ 4 M yr~'. At r 0.1 x 10 16 cm, i.e. at 
the inner boundary of the envelope, density and velocity are set 
to the actual stellar wind properties for age zero, M = 2.5 x 
10~ 7 M© yr 1 and v = 60 km s _I (for Z = Z GD , see text). 



content below 1/10 solar ( Pequignot & Tsamis|2005| Sandin et 
|al.|2010||Stasiriska et al.|2010) . 

The wind is treated in the same manner as originally devised 
by |Marten & Schonberner| ( |1991| l, but its dependence on metal- 
licity, more correctly on the elements C, N, and O, is approxi- 
mately accounted for by simple c orrection factors, i.e. M oc Z 69 
dVink etal.|200ll and v^ocZ 013 dLeitherer et all 1992k. The out- 



come is a nearly linear dependence of the wind luminosity with 
metallicity, L wind = 0.5 M i& oc Z 95 . 

The surface of the star is always assumed to radiate like a 
black body, i.e. a dependence of the radiation field on the chem- 
ical composition of the stellar atmosphere is, for simplicity, not 
considered here. The influence of the metal content on details of 
the stellar evolution is expected to be small and is thus also ne- 
glected. Throughout this pilot study we employed for all metal- 
licity cases our standard set of central-star models already used 
in our previous work. For the wind power calculations, the cen- 
tral star is assumed to have the same metal content as the nebular 
gas. 

In every case we used the same initial envelope structures, 
independent of their metal content, thereby ignoring the pos- 
sible influence of the metallicity on final AGB mass-loss rates 
and wind velocities. Our present observational and theoretical 
knowledge about the AGB wind properties as a function of 
metallicity is still too meager as to allow a more detailed de- 
scription of initial configurations in terms of their metal content, 
although progress, theoretical and observational, is encouraging 
([Marshall et al.|2004||Mattson et al |2008l|Sandin|2008l[S"a ndin 
l et al.p010| |van Loon et al.||20051 |Wachter et al.||2008| |Wood 
|etaL][1992> . There are clear indications that the wind speed is 
somewhat reduced at lower metallicity, but the dependence of 
the mass-loss rate on metallicity, and on stellar parameters, re- 
mains still unclear. It must also be noted that the particular evo- 
lutionary stage during which the remnant leaves the AGB, which 
is also the only relevant here, is not covered by any of these stud- 
ies. 

For a first assessment of the influence of the metal con- 
tent, we selected a circumstellar envelope structure with initially 
a = 3 and v — 10 kms -1 from our set of initial models with 



power-law radial density profiles, p oc r a , used in Paper II 
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Fig. 3. Evolutionary path of the 0.595 Mq post-AGB model used 
for the a — 3 hydrodynamical sequences of this work. Post-AGB 
time marks (circles) are separated by 10 3 years. 



There it was shown that the expansion speeds seen in 
Galactic disk PNe with round/elliptical shapes imply radial den- 
sity gradients with power-law exponents between 2.5 and >3, 
where the larger exponent is typical for the more evolved objects. 
The choice of a density gradient with a = 3 for all metallicities 
reflects just a compromise in order to make the model grid as 
simple as possible without losing its significance for practical 
applications. 

Figure[2]illustrates the structure of the initial models in some 
detail. This envelope model was coupled to a hydrogen-burning 
post-AGB model of 0.595 M0 used in Paper II whose evolution 
in the Hertzsprung-Russell diagram is depicted in Fig. [3] The 
actual variations of the central star wind, based on the adopted 
metallicity range and the approximations introduced above, is 
illustrated in Fig. [4] There are two facts which emerge from this 
figure: 

(i) The wind power is very weak at the early phases of evolution 
but increases by about two orders -of-magnitude (L w j n d 00 
before it decreases again with the fading central star. 

(ii) Even for the case of the most powerful wind (3 Zgd) its me- 
chanical power is only a very small fraction of the photon 
luminosity. 

The evolution of the whole system, star and circumstel- 
lar envelope, is followed across the Hertzsprung-Russell dia- 
gram towards the white-dwarf cooling path, employing our ID 
radiation-hydrodynamics code. An important aspect of our hy- 
drodynamical simulations is that we can switch off the hydrody- 
namics by setting the gas velocity to zero at any time in the sim- 
ulation. The model is then able to settle into its equilibrium state 
for fixed density profile and radiation field which is reached if 
the electron temperature distribution remains constant through- 
out the entire computational domain. This kind of model corre- 
sponds then to a standard photoionisation model, but with the 
additional feature of having a velocity field and also a halo. We 
are thus able to estimate in a self-consistent manner possible 
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Fig. 4. Wind luminosities, L w i n d = 0.5 M v^, for different metal- 
licities as a function of the (post-AGB) age, compared with the 
photon luminosity of the 0.595 M0 post-AGB model shown in 
Kg.|3] 

Table 2. Overview of the model sequences used in this work. 
The table lists the central star masses (Col. 1), the AGB-wind 
properties (Cols. 2-A), the different metallicities covered for 
each model (Col. 5), and references to detailed discussions of 
the inital models (Col. 6). 
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Type 


Magb 




Metallicity a 


Reference 


[M©] 




[Moyr 1 ] 


[km s" 1 ] 


[Z GD ] 






0.595 


a = 3 


1.3 x 10- 4b 


10 


3... 1/100 


Paper II 


0.595 


Hydrodynamics sim. 


^12 
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Paper 111 


0.605 


Hydrodynamics sim. 


^12 


3,1 


Paper IV 


0.605 


a = 2 


1.0 x 10~ 4 
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Paper 1 
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3. ..1/10 


Paper 1 
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a = 2 


0.5 x 10~ 4 
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1/3 






0.696 
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1.0 x 10~ 4 
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Paper I 




0.696 


a = 2 


1.0 x 10~ 4 


20 
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a = 2 


1.0 x 10~ 4 
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0.696 


a = 2 


0.5 x 10~ 4 


15 


1 






0.696 


a = 2 


2.0 x 10~ 4 


15 


1,1/3 







a Metallicity ranges are are filled in steps of 0.5 dex. 

b Value at a distance r = 1.0 x 10 16 cm from the central star. 



systematic errors introduced by imposing equilibrium conditions 
for systems which are in reality influenced by dynamics. 

Central stars of about 0.6 M0 are not luminous enough as 
to harbour the brightest (and thus easiest observable) planetary 
nebulae in a distant stellar system. We supplemented thus the 
already existing hydrodynamical sequences with central stars of 
0.625 and 0.696 Mq, both with Z = Zgd, which were discussed 
in [Paper T] and Schonberner et al. (2007 Paper IV hereafter), 
by additional sequences with other metal contents: Z = 6 Zgd 
(0.696 M© only), 3Z GD , Z GD /3, and Z GD /10. Except for the 
metallic ity, the initial models are identical with those used in 

71 „. _ o ...;.u ji> _ iri-4 iv/r_ „ . _ it 1™, „-l 



Paper ha -2 with M agb = 10~ 4 M yr , v. dgh = 15 kms" 1 . In 
these cases the choice of a constant mass outflow (i.e. a — 2) 
appears quite reasonable since these rather massive central stars 
have a very short post-AGB phase during which the PN remains 
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Fig. 5. Density, velocity and temper- 
ature profiles of models from the 
a = 3 sequences, selected after ^3500 
years of post-AGB evolution, and for 
the metallicities indicated above each 
panel. The central star is at the origin, 
and its parameters are T c g =: 41 000 K 
and L = 5534 Lq. The velocity oscil- 
lations seen in some early models are 
due to a temporary numerical instabil- 
ity in the hot bubble caused by our 
heat conduction routine, but have no 
impact on the model dynamics. Upper 
row: the radial profiles of heavy parti- 
cle densities (thick solid), electron den- 
sities (dotted), and gas velocities (thin 
solid). The particle densities are nor- 
malised and must be multiplied by the 
factors given in the individual panels 
to get the true densities. Lower row: 
the radial run of the electron tempera- 
tures, where the dotted lines represent 
the equilibrium temperatures (see text 
for the details, and also Sects. [3~L2] and 
[3T3J. 



very compact and is expected not to 'see' larger variations of the 
circumstellar density gradient while it is expanding. 

In order to estimate possible influences of different, albeit 
constant, mass-loss rates and outflow velocities on the nebular 
kinematics we computed several additional sequences for the 
0.696 Mq post-AGB model with AGB mass-loss rates and wind 
velocities changed. For the 0.625 M0 post-AGB model, we com- 
puted one additional Zgd/3 sequence with M ag b = 0.5 x 10~ 4 
M0 yr _1 and u ag b = 7.5 km s _1 . An overview of all the sequences 
used in this paper is presented in Table [2] 

In the following sections we discuss in more detail implica- 
tions which follow from our simulations and which are impor- 
tant in interpreting PNe in distant stellar systems with varying 
metallicities. We will thereby distinguish between the a - 2 and 
a — 3 model sequences. 

3. Results 

3.1. The a = 3 sequences 



3.1 .1 . Structures and kinematics of the models 



Figures. [5[j9] illustrate how the general evolution of the model 
nebulae is influenced by their metal content. All sequences have 
the same central star of 0.595 Mq and start with the same initial 
model depicted in Fig. [2] The figures show 5x4 snapshots along 
the stellar track, viz. for five positions and for four metallicities 
each, Z = 3Zgd, Zgd, Zgd/10, Zgd/100. The positions selected 
are at the early ionisation phase at about T e g - 40 000 K (Fig. [5}, 
at the high-ionisation stage at about 100000 K (Fig. [6]), at the 
maximum stellar effective temperature (Fig. [7J, during nebular 
recombination at a stellar luminosity of about 500 Lq (Fig. |8}, 
and at an even lower luminosity of only 240 Lq after reionisation 
has started (Fig. [9]). Note the different scaling of each individual 
panel within each figure, and also the different radial ranges of 
the figures. 

An extensive discussion of the principles of PN evolution 
based on radiation-hydrodynamics simulations has been given 
in Paper I, which the reader is referred to for details. We concen- 



trate here only on those aspects relevant in connection with the 
metal content of the models. 



The upper panels of Figs. [5jj9] demonstrate how structure and 
kinematics evolve with time. We see a strong trend with metal- 
licity: the lower the metal content, the larger and thicker be- 
come the nebular shells and the smaller the wind-blown cav- 
ities. Responsible for this behaviour are two factors, both de- 
pending on the metal content. The first is the expansion speed of 
the hot, ionised gas which scales with the sound velocity oc yT^ 
(cf . [Paper II| l and becomes higher at lower metallicities because 
of the reduced line cooling efficiency. The second is the wind 
power which decreases with metallicity as outlined in Sect. [2] 
(Fig. |4| and loses its ability to compress and accelerate the in- 
ner nebular parts. This is seen in the two metal -poor sequences 
where the nebular density falls off radially more gradually with a 
nearly linear slope (Figs.[5}{8] top panels). A sharp density con- 
trast between the rim and the shell, typical for the metal -richer 
models, is therefore an indication of a strong wind. 

Once the central star has passed the position of its maxi- 
mum effective temperature in the Hertzsprung-Russell diagram 
(Fig. [7J and starts to fade quickly, also the wind power becomes 
progressively weaker in all cases (see Fig.|4]i. Figure|8]displays a 
moment during the recombination phase as the central star fades, 
and Fig. [9] a stage after reionisation has begun to dominate the 
inner nebular regions, forcing there the matter to develop a posi- 
tive density gradient, very similar to the typical situation behind 
a D-type shock during the first ionisation. 



Structures Figure [T0| summarises the complete size evolution 
with time for all six metallicities investigated and illustrates 
clearly the importance of the wind interaction, given by the pres- 
sure of shocked wind (bubble) gas, as compared to the thermal 
pressure of nebular matter, and how these pressures change dur- 
ing the course of evolution. From this figure it becomes evident 
that wind interaction is not the main driver for the nebular ex- 
pansion, although its importance becomes larger during the end 
phases of evolution and generally with metallicity. Rather, the 
expansion is initiated by thermal pressure differences, caused by 
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Fig. 6. The same as in Fig. [5] but after 
-7150 years, the 3Z GD model is now 
also optically thin. The stellar param- 
eters are T cff 100900 K and L = 
4580 L©. 
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Fig. 7. The same as in Fig. [5] but af- 
ter ^9840 years. The stellar parameters 
are now r eff = 146 870 K and L = 
1845 Lq, corresponding to the turn- 
around point of the stellar track seen in 

Fig.0 



photo-heating of the former neutral circumstellar gas, and con- 
tinues independently of the stellar wind. Wind interaction is only 
responsible for compressing the inner parts of the expanding gas 
and preventing the shell from collapsing. This last statement is 
also true during the final evolution along the white-dwarf se- 
quence: Even at the lowest metallicity considered here the hot 
bubble is still expanding, although the stellar wind power is re- 
duced by about one order-of-magnitude (cf. Fig. [4]). 

For the metal-richer models with low expansion rates and 
more powerful central-star winds, wind interaction becomes 
more important and leads to very (geometrically) thin nebu- 
lar shells. However, even for the highest metallicity considered 
here (3Zgd) the double-shell structure remains visible during the 



whole high-luminosity part of evolution. The rim succeeds in 
overtaking the shell not before the end of evolution (cf. Figs. [7] 
to [9). Neglecting photo-heating would lead to single-shell (i.e. 
rim only) structures in all cases. 

From the observer's point of view, only the relative sizes 6R 
(also shown in Fig. 10 1 are of interest because only these can be 
measured distance-independently. In general, the rapid increase 
of 5R at the beginning of the evolution reflects the dynamics of 
the increasing thermal pressure due to ionisation, and the de- 
crease of SR later on is due to the increased nebular size and the 
stronger stellar wind. At the lowest metallicity, Z = Zgd/100, 
the relative sizes are close to 0.9 for most of the time, in con- 
trast to the metal-rich models, Z = 3Zqd, which are much more 
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Fig. 9. The same as in Fig. [5] but af- 
ter ^13 850 years. The stellar param- 
eters are now T e g ^ 119 800 K and 
L = 240 Lq, typical for the reionisation 
stage when the nebula expands around 
a central star with virtually constant lu- 
minosity and temperature. 



compressed and whose relative thicknesses do not exceed «0.6. 
The Galactic disk composition favours medium thick objects 
with 6R between ^0.5. . . 0.7, consistent with the observations of 



round/elliptical PNe (cf. Fig. 6 in Paper IV I. 

The very different density profiles which our models develop 
during their course of evolution are, of course, also reflected in 
their surface brightness distributions. The case is illustrated in 
Fig. 



11 



where the Hp surface brightness profiles of the models 
shown in Figs. [5}|7]and[9] are displayed. 

One sees a clear trend with metallicity: The metal-rich mod- 
els (Z S Zgd) develop a deep central cavity and a pronounced 
rim-shell structure where the shell brightness is only about 10 % 
or less of the maximum rim brightness. The metal-poor models 



display a more gradual brightness decline with distance from the 
star, with a nearly linear (negative) slope for the lowest metallic- 
ity investigated (Z = Zgd/ 100). The rim-shell dichotomy known 
from most PNe is no longer existent. Additionally, the central 
cavity becomes much smaller and is nearly filled up by the pro- 
jection effect. 

Kinematics The velocity field is very similar in all cases. Once 
ionisation has started, the gas velocity increases nearly linearly 
with radius. The post-shock speed (i.e. the gas velocity imme- 
diately behind the outer edge of the nebular shell) increases 
with time in line with the shock acceleration (Figs. |5j|9|. Note 
that during the early stage of nebular evolution the innermost 
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Fig. 10. Development of nebular sizes with time for the a = 3 
sequences with various metallicities. Right ordinate: inner (R c d) 
and outer radii (R ou t), with the different metallicities indicated 
by various (overlapping) gray shades; see legend. Left ordinate: 
relative thicknesses, 6R = (R ou t - R c d)/Rout, f° r the same se- 
quences. Metallicity decreases monotonically from the bottom 
curve towards the top curve. 



part of the ionised shell expands slower than the former AGB 
wind because it is decelerated by the high thermal pressure. This 
fact holds also during the whole evolution in the low-metallicity 
cases (top panels in Figs. [5]-[9]l. The post-shock speed is usu- 
ally also the maximum expansion velocity within the whole neb- 
ula, with the exception of the more metal-rich models with their 
strong stellar winds which create dense rims and accelerates 
them to velocities greater than the post-shock one (Figs. [7Jj9]i. 
An extreme case occurs for Z = 3Zgd where the rim becomes 
very dense, thin, and fast, and swallows finally the outer shell 
(see Figs. |7}{9). 

The extremely metal-poor models reach quite high post- 
shock (gas) velocities, up to 60 kms~', which is twice the value 
achieved by the models of the 3Zgd sequence. However, the 
most conspicuous features seen in Figs.[3J|9](bottom panels) are 
the strong shocks that these models develop with time where the 
shell interacts with the largely undisturbed AGB wind. The post- 
shock temperatures reach extremely high values, up to 50 000- 
60 000 K, and the whole post-shock thickness can be substantial 
for the most metal-poor cases. Despite of this no observed sig- 
natures are produced by these shocks because of the very low 
matter densities involved. 

The overall nebular expansion properties and their depen- 
dence on metallicity is summarised in Fig. [12] (left) where the 
propagation speeds, R out , of the shock, which defines the outer 
edge of the model PNe, are plotted. Starting at rather low values 
the shock propagation speed increases in all cases very rapidly 
until the nebular shell becomes optically thin. Afterwards the 
shock speed (relative to the upstream flow) is determined by the 
density profile, p oc r~ 3 , and the electron temperature (or sound 
speed) which becomes higher as the central star gets hotter and 
the ionising photon flux more energetic (cf. Paper II Sect. 3.1 
and Fig. 7 therein). For instance, the sudden velocity increase 
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Fig. 11. Normalised radial surface brightness (intensity) profiles 
in H/3 of the models shown in Figs. [5|j7]and[9] For clarity, the 
inset in the top left panel shows the surface brightnesses at an 
enlarged scale. 



between 6000 and 7000 years seen in all sequences is due to the 
second ionisation of helium, causing a fast growth of the electron 
temperatures (see also Fig. 18 1. 

As expected, the shock velocities increase with decreasing 
metallicity, reflecting the values of the nebular electron tempera- 
ture (cf. Figs.[5]-[9] and also Fig.[T8|. The effect becomes smaller 
for the lowest metallicities just because also the electron tem- 
perature increase levels off somewhat (cf. Fig. 18 below). The 
total range of R out at the end of our simulations (^15 000 yr) 
is ^55-80 kms~' for the metallicities used here. In the 3Zgd 



sequence with its powerful stellar wind we see the rim shock 
overtaking the more slowly expanding (38 kms ) outer shock. 
This occurs at t 13 000 yr when Rr^, jumps from 38 km s _I to 
^out = ^iim — 60 km s _1 (cf. Figs.[7ll9l leftmost panels). 

The maxima of the gas velocities are always achieved right 
behind the outer shock, and are depicted in the right panel of 
Fig. 12 The trend with metallicity is, of course, the same as 
seen for R out , but the absolute values are lower: starting from the 
shared AGB wind velocity of 10 km s , the gas becomes rapidly 
accelerated to maximum values of 32-62 kms -1 , depending on 
metallicity. We see again the velocity jump in the 3Zgd sequence 
when the shell becomes swallowed by the faster rim (cf. Figs.[7|- 
[9] leftmost panels). 

Our simulations demonstrate clearly that the expansion of 
a planetary nebula is ruled by the electron temperature of the 
shell gas (cf. Paper II), and not by the wind from the central 
star as is predicted by the favourite theory of interacting winds 
put forward by |Kwok et al.| (fl978). Wind interaction is only re- 
sponsible for the shape and acceleration of the rim, i.e. the inner, 
bright parts of a PN, and under most conditions of our simula- 
tions the rim expands slower than the shell. With the reasonable 
assumption made here that the stellar wind strength decreases 
with metallicity, it follows that the nebular expansion of the outer 
shell is fastest if the wind is weakest! 

The question arises whether the large variation of the expan- 
sion velocity with nebular age and metallicity can be measured. 
For spatially resolved objects there is the possibility to mea- 
sure the fast moving matter behind the outer shock front using 
high-resolution line profiles taken for the central line-of-sight 
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Fig. 12. Expansion properties of the 
model nebulae of the six a = 3 se- 
quences with different metallicities vs. 
the post-AGB age. The five black dots 
along the abscissas refer to the snap- 
shots depicted in Figs. |5j|9] The initial 
velocities of the models are 10 kms~' 
(cf. Fig.[2]l. Left: propagation rates of 
the outer shocks, R oul , which define 
the nebular sizes. The velocity jump of 
the 3Zqd model at t = 13 000 yr marks 
the moment when the faster expanding 
rim starts overtaking the outer shock, 
and R out = R Ilm afterwards. Right: the 
gas velocities immediately behind the 
outer shock. Again, the velocity jump 
of the 3Zgd model, now already at 
t m 12 000 yr, is due to the disappearing 
shell. 



Fig. 13. Expansion properties as de- 
rived from the half width at half max- 
imum (HWHM) of the spatially in- 
tegrated line profiles for the same six 
a = 3 sequences shown in Fig.[l2] vs. 
the post-AGB age, for H/3 (left) and 
[Om] 5007 A (right). Again, the five 
black dots along the abscissas refer to 
the snapshots depicted in Figs. [5}j9] 
The shadowed areas indicate the ther- 
mal Doppler broadening, computed by 
means of the mean electron temper- 
atures (see next section). The lower 
boundary refers to the models with 
the highest metallicity, 3Zqd, the up- 
per boundary to those with the lowest 
metallicity, Zqd/100. 



(Corradi et al. 2007). In the case of distant objects which cannot 
be spatially resolved, only the half width of the (integrated) line 
profile can be used. We computed therefore spatially integrated 
line profiles for H/3 and [O in] 5007 A and determined Vhwhm as 
a measure for the nebular expansion. The results are summarised 



in Fig. 13 



First of all, the spread in Vhwhm due to the metal content is, 
if compared with the situation in Fig. 12 surprisingly small and 
partly irregular, and the increase of Vhwhm with time appears to 
be more gentle in both lines, especially during the early phase of 
evolution. This can be understood as a consequence of integrat- 
ing over the whole object which gives much more weight to gas 
elements with low velocity components along the line-of-sight 
if compared to the central line-of-sight case. This effect is en- 
hanced by the density structure of the models: the gas velocity is 
lowest where the gas density is highest. 

The HWHM velocities based on H/? are higher than those 
derived from 5007 A during the early evolution, which is due to 
both the different thermal line widths, indicated in Fig.[l3]by the 
shadowed areas, and the ionisation structure: if the central star 
is still not very hot, the 2+ zone is confined to the inner part 
of the ionised region only, which is strongly decelerated by the 
thermal pressure (cf. Fig. 151. Thus, Vhwhm as measured from 
the (integrated) [O in] line starts at very low values, which are 
significantly below the original AGB-wind velocity of 10 km s _1 
assumed here. Later on, at larger ages when the Doppler width 
of HjS becomes larger than the thermal one and the 2+ zone 



is extending to the outer nebular edge (i.e. to the outer shock), 
both lines behave very similarly. In the following all velocities 
are from [Om] only if not specified otherwise. 

In any case, the HWHM method severely underestimates the 
true expansion of the nebula which is given by the propagation 
of the outer shock front, as seen in Fig. [12] (left panel). For in- 
stance, at age 5000 years the HWHM velocities are between 10 
and 20 km s _1 while the true expansion velocities are already sig- 
nificantly higher, viz. between 30 and 60 km s _1 ! Thus, the rather 
large sensitivity of the nebular expansion on metallicity is almost 
lost by using spatially integrated line profiles. An observational 
test of the predicted dependence of the expansion velocities on 
metallicity by using an appropriate sample of PNe drawn from 
galaxies with different metal contents as done by Richer| ( 2006) 
appears to be difficult if not impossible. We will come back to 
this point in more detail in Sect. |4] 



Applications It is tempting to apply our model sequences to 
metal-poor objects of the Milky way. The two well-known ob- 
jects NGC 4361 and NGC 1360 appear especially suited for this 
purpose because structure and velocity informations are avail- 
able from the literature. NGC 4361 belongs to the Galactic halo, 
while NGC 1360 is not known as an halo object but is metal- 
poor as well (see below). According to Mendez et al. (1992), 
both objects consist of a very highly excited nebula surrounding 
a very hot and luminous central star. 
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Fig. 14. Expansion properties as de- 
rived from the line peak separation 
of spatially resolved line profiles for 
the same six a = 3 sequences shown 
in Figs.[l2] and \H\ vs. the post-AGB 
age, for H/3 (left) and [Ora] 5007 A 
(right). The line profiles are simulated 
using a central numerical aperture of 
1 x 10 16 cm (or Cf.'67 at a distance of 
1 kpc). At low ages, the profiles are 
still singly peaked, and thus V pcak = 0. 
The [Om] profiles are only plotted up 
to model ages of about 10000 yr be- 
cause beyond this age the profiles de- 
velop a complex structure due to recom- 
bination. The five black dots along the 
abscissas correspond again to the snap- 
shots depicted in Figs.[5}j9] 
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Fig. 15. Line profiles of Hf3 (top) and [Om] 5007 A (bottom) 
computed for the models displayed in Fig. [6] for a central aper- 
ture of 1 x 10 16 cm and broadened by a Gaussian of 10 kms -1 
FWHM. The profiles are normalised to / = 0.25 and shifted 
by A/ = 0.25 each for more clarity. Different symbols mark 
shock (open circles) and post-shock velocities (filled circles) of 
the models. 



Since both objects are spatially resolved, we computed the 
line profiles for an aperture centred on the position of the central 
star (Fig. [14} . At some time the lines become double-peaked, and 
this figure demonstrates that the expansion rates as derived from 
the line peak separations, y pea k, are very similar to those derived 
from the integrated profiles: The dependence on metallicity is 
weak or even partially absent (hydrogen lines only), and the true 
expansion velocity is strongly underestimated as well, especially 
for the metal-poor cases. 

For illustration, the (resolved) line profiles of hydrogen and 
[Om] for model ages of about 7100 years (cf. Fig. [6]i are ren- 
dered in detail in Fig. 15 At these ages, our models have very 
high nebular excitations and still luminous central stars with 



about 100000 K effective temperature. The peak line emission 
corresponds to the denser inner regions where the emission mea- 
sure is high but the gas velocity still rather modest (cf. upper 
panels of Fig. |6j. The fast moving nebular layers immediately 
behind the outer shock contribute very little to the total profile 
and are masked by the strong emission from the denser nebular 
regions. This figure emphasises vividly how the real expansion 
of a PN is underestimated by employing the peak separation of 
resolved line profiles]^] 

Figure [l5]displays also an interesting line width behaviour: 
The widths increase with decreasing metallicity. This is due to 
the different ranges of gas velocities encountered in the models 
(cf. Fig. |6[). From the highest to the lowest metal content shown 
in Fig. [6Tthis velocity range varies from about 10 to about 44 
kms -1 . Also, we see for the lower metallicities a trend that the 
[O in] lines have larger peak separations than the hydrogen lines. 
The reason is the ionisation profile of oxygen: At the evolution- 
ary stage considered here, the 2+ concentration increases to- 
wards the outer nebular region where the velocities are higher 
(cf. 6th left panel in Fig. [20 1. 

Guided by our simulations, we expect metal-poor PNe to 
show a more smooth or diffuse morphology without a clear-cut 
signature of a central cavity (or hole) and a sharp outer bound- 
ary. The surface brightness profile of NGC 4361 shows indeed 
neither a pronounced central cavity nor a bright, distinct rim 
( |Monreal-Ibero et al.|2005|l. The metal content of this object is, 
on the average, O.IZgd ( Howard et al.|1997 Table 4 therein). 
It has also one of the highest electron temperatures, 19 300 K, of 
the whole [Howard et aT] sample (see Fig. [1). The central star has 
an effective temperature of 82 000 K (M endez et al.|1992 1. 

In Fig. 16 we compare the Ha surface brightness of NGC 
4361 with the predictions of our a — 3 models discussed here. 
The surface brightness is taken from a cut along the semi-minor 
axis in order to avoid complications due to the weak ansae seen 
at the poles. We conclude from the figure that only the models 
with a reduced metal content, Z Zgd/10, provide a quite rea- 
sonable agreement with the observation in terms of central dip 
and intensity profile. 

The kinematics of our models can be tested as well: 
The latest velocity measur ements of NGC 4361 are 



V P eak(5007) = 27 ± 2 kms" 1 ( |Muthu & Anandarao 



2001 



2 For the 3Z G d case the post-shock velocity coincides just by chance 
with the peak velocity because the rim which is resposible for the peak 
emission is being accelerated by the strong wind/bubble pressure (see 
upper left-most panels of Figs. [6] to [5} . 
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Fig. 16. Ha surface brightness of NGC 4361 (dots, left panel) 
derived from a cut along the semi-minor axis of the Ha image as 
indicated (right panel). We extracted this image from observa- 
tions reported in Monreal-Ibero et al. (2005 ). The dark horizon- 
tal bars are artefacts of the VIMOS spectrograph. Also shown 
are the predicted Ha surface brightness distributions of models 
with different metal contents (see figure legend), taken from the 
a = 3 sequences. These models are chosen such that they match 
the observed size of NGC 4361, assuming a distance of lkpc. 
The model ages vary from about 6700 years (Zqd/30, fast ex- 
pansion) to 9100 years (Zgd, slow expansion). 



and V peak (Hff) = 22 + 5 kms' 1 , V pea k(5007) = 26 + 5 kms" 1 
( |Medina et al.||2006) . |Mufhu & Anandarao] ( |2001| l provide 
[O m] line profiles, and the one taken from the central region of 
NGC 4361 (their Fig. 4c) is compared with theoretical profiles 
generated from the models of Fig. 



16 Our metal -poorer models 



predict the right observed apparent expansion velocity, broader 
line profiles, and also V pe ak(5007) > V pea k(Hy8). 

The other case, NGC 1360, cons ists also of a hot (97 00 K), 
luminous (^5000 L0) central star (Traulsen et al. 2005) and a 
rather smooth, extended nebula without a clear cavity/rim/shell 
signature (Goldman et al. 2004, Fig. 5 therein). It is reported 



Sect. 



to have a very high electron temperature (16 500 K, Kaler et al. 
1990), indicative of a correspondingly low metal content (see 
3.1.2| i. Abundance determinations for thi s object are rare: 



Kaler et al. provide only e(0 ) = 7.8, while |Manchado et al. 



( |1989| l found e(O) > 8.2, e(Ne) > 7.6, e(Ar) = 5.8, which are, 
however, based on an electron temperature of only 10000 K. 
A higher electron temperature as mentioned above would re- 
sult in correspondingly lower abundances. Traul serTet al.| (f2005 ) 
provide also photospheric abundances of the central star which, 
converted to the notation used here, are: e(C) = 8.3, e(N) = 7.7, 
and e(O) = 8.3. Thus, NGC 1360 is certainly also a metal-poor 
object (compare with the Galactic disk abundances in Table[T| 



Garcia-Diaz et al.|p008) measured recently the apparent ex- 



pansion of NGC 1360 by means of high-resolution echelle spec- 
trograms: V pe ak(Hor) = 26 km s _1 , consistent with older measure- 
ments of 24 kms~' by |Goldman et al.| ( |2004| . Again, the agree- 
ment with the predictions of our models for this particular stage 
of evolution (^7100 yr, Fig. 14 1 is satisfying. We repeat that the 
real expansion rate of both objects is expected to be considerably 
higher, viz. up to a factor of about two (cf. Fig. 12 left panel, and 



Goldman et al. (2004i determined also an elliptical spatio- 
kinematic model of the nebula structure of NGC 1360, based 
on the position-velocity ellipse and the surface-brightness distri- 
bution. The basic result is that the gas density falls off outwards, 
while the gas velocity increases outwards, i.e. the density (or sur- 
face brightness) is highest close to the central star where the ex- 
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Fig. 17. Observed central 5007 A line profile of NGC 4361 
(connected crosses) as measured by M uthu & Anandarao] (2001) 
compared with the profiles generated from the models used in 
Fig. [16] As in Fig. 15 the theoretical profiles are computed for 
a central aperture of of 1 x 10 16 cm and then broadened by a 
Gaussian of, now, 12 kms" 1 FWHM. 



pansion velocity is lowest. The fastest expanding shell has ^35 
kms 1 at the equator and ^70 kms~' at the pole. The density, 
velocity, and surface brightness profiles predicted by our metal- 
poorer models (Fig. [6] top panel s, and Fig. [TT|l a re in agreement 
with the mean properties of the [Goldman et al. model of NGC 
1360._ 

constructed a dynamical (2D) 



Garcf a-Dfaz et al.| (2008 



model including magnetic fields, but their main point was the 
following: They assumed that the stellar wind has stopped 
shortly after the formation of the PN (i.e. after 1000 years) and 
followed the dynamical evolution with a collapsing hot bubble 
for additional 10000 years. Despite the authors' claim that their 
model "successfully reproduce many of the key features of NGC 
1360", we note following inconsistencies: 



1. 



The final nebular model at t 
p008 



Diaz et al. 



1 1 000 yr in Fig. 6 of |Garcia- 
shows a positive gradient for both the ve- 



locity and density: the outer fastest layers are also the dens- 
est, i.e. the model is limb-brightened (except at the poles), 
which is in clear contradiction with the observations. 
Despite the long simulation time of 1 1 000 years, an evolu- 
tion of the star is obviously not considered. The stellar model 



envisaged is rather massive (^0.75 M0, Vassiliadis & Wood 
1994) and crosses the HR diagram within only about 1000 
years and fades thereafter, accompanied by a drop of wind 
power (cf. Figs. [3] and [4j. After 1 1 000 years, the model's 
luminosity is then only a; 100 Lq. But presently the central 
star of NGC 1360 is still a high-luminosity object with sev- 
eral 1000 L.0, luminous enough to sustain a radiation-driven 
wind which, because of the low metallicity, may be too weak 
for developing spectroscopic signatures. 

We emphasise at the end of this discussion that models taken 
from the a — 3 sequences displayed in Figs. [5] to [9] can explain 
the (mean) properties of two rather well-known metal-poor ob- 
jects very successfully, although they were not designed to fit 
any particular object! Furthermore, according to our models, it 
is not justified to conclude, just from the non-existence of an ap- 
parent central hole in the surface brightness distribution of a PN, 
that the central-star wind has already stopped and that the nebu- 
lar gas is now falling back onto the star. Rather, a detailed analy- 
sis of the density and velocity structure of the object in question 
appears to be mandatory in order to decide whether a deceased 
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central-star wind is responsible for a more smooth or diffuse ap- 
pearance of a PN, or whether we observe a low-metallicity sys- 
tem with a weak central-star wind and a more extended nebular 
structure. 

We acknowledge, however, that our spherical models are cer- 
tainly too simple in order to provide detailed models of the ob- 
jects discussed in this section. Important further ingredients are, 
e.g., inhomogeneities, jets, and other means that impose depar- 
tures from sphericity. We note also that magnetic fields have 
been detected in NGC 1360 by |Jordan et al.| ( |2005) , but whether 
they are really important in shaping NGC 1360, as believed by 
|Garcia-Diaz et al.| (|2008 ), must be seen in the future when more 
realistic simulations become possible. 

3.1.2. Electron temperatures 

We computed the mean electron temperatures of the models ac- 
cording to 



(T e ) = 



J T e (r)N e (r)Nj(r)dV 
jN e (r)Ni(r)dV 



(1) 



with N e (r) being the electron number density, N[(r) the +2 
number density, T e {r) the electron temperature, and dV the vol- 
ume element. The upper boundary for the volume integration is 
set at the outer shock front; the halo (= ionised former AGB 
wind) is thus excluded. In passing we remark that Eq. ([T]i corre- 
sponds in practice to an emission weighted mean of the electron 
temperature, in contrast to the usual meaning of a mean value. 
Differences are only expected for cases with radial temperature 
gradients and/or strong shocks as seen in the late, metal-poor 
models of Figs. [3}j9] 

The determination of the mean electron temperature from the 
model structure by using Eq. ([TJ is convenient but different from 
the one used in practice if the nebular structure is not known. 
The alternative is to compute first the (total) line strengths from 
the model and to apply then, e.g., Eq. (5.4) of Osterbrock ( 1989 ) 
to get a volume averaged 7om from the ratio Rom. Provided 
the collision strengths are correct, both methods should give 
the same result. We performed a test using two sequences with 
rather extreme metallicities, viz. with Zgd and Zgd/100. In the 
first case we found always \(T e ) - 7b ml ^ 60 K, in the second 
|(T e ) - 7om| ^ 170 K. Thus the (relative) differences between 
both methods are less than about 1%, and the temperature de- 
termination according to Eq. ([T| provides an excellent value for 
?o in- 

The run of the mean electron temperature, (7 e ), is displayed 
for all a — 3 sequences in Fig. 18 and we see, as expected, a 



strong temperature dependence on metallicity. The first rapid 
increase of the mean electron temperature is, of course, due to 
hydrogen ionisation. Already at this stage one sees clearly how 
rapidly the electron temperature increases with decreasing cool- 
ing efficiency. Then the temperatures increase further as the cen- 
tral star becomes hotter and its photons more energetic. The sec- 
ond temperature "jump" after about 5000 years is due the sec- 
ond ionisation of helium, providing additional energy input to 
the gas. 

Because the effective temperature decreases once the cen- 
tral star begins to fade rapidly at an age of about 10000 years 
(cf. Fig. B), the electron temperature in the nebula must drop 
accordingly. This is clearly seen in Fig. [18] (left) for the more 
metal-rich models in which expansion cooling is not very im- 
portant. However, the relation (7 e )-7 e ff differs between the stel- 
lar high- and low-luminosity branch because of differences in 



the nebular ionisation: At a given stellar temperature, the neb- 
ular ionisation is higher and line cooling lower at the high- 
luminosity branch. Hence, the electron temperature along the 
low-luminosity branch of evolution is significantly below that 
along the corresponding high-luminosity branch, independently 
of non-equilibrium effects. 

An interesting effect is seen for the more metal-poor se- 
quences (Z < Zgd/3) which expand faster and become hence 
more dilute: Here expansion cooling becomes relevant during 
more advanced evolutionary stages and forces the nebular elec- 
tron temperature to drop at progressively lower stellar tempera- 
tures (cf. Fig. 18 right panel). A more detailed discussion about 
expansion cooling is given in Sect. |3.1.4| 

Note that the interplay between line and expansion cooling, 
and their relative importance, depend sensitively on gas density 
and expansion velocity. It is, however, clear from our simula- 
tions that a low metal content results in high electron tempera- 
tures because of less line cooling, leading consequently to faster 
expanding and more diluted nebulae which are then prone to de- 
partures from thermal equilibrium. The grey region seen in the 
right panel of Fig. 18 delineates the approximate area occupied 
by (r e ) of all the hydrodynamical sequences previously used in 
Paper I and Paper IV which have very different combinations of 



initial AGB-envelope configuration and central star, but all with 
the Zgd composition. These models sequences demonstrate that 
only a low metal content can be responsible for unusually high 
electron temperatures in nebular envelopes. 

We note also that the conspicuous shock regions seen in Figs. 
[5}|9]do not contribute to the electron temperatures, although the 
post-shock temperatures become very high and are included in 
the integration according to Eq. (QJ. The reason are the low 
ion (i.e. 2+ ) densities, giving these outer shock regions only 
a minute weight in the integral. 

3.1.3. Comparisons with equilibrium models 

With the computation of equilibrium models (with respect to 
ionisation and thermal energy) at selected positions along the 
evolutionary path we have a unique tool to investigate under 
which circumstances non-equilibrium effects may become im- 
portant. 

Electron temperatures In equilibrium, the electron tempera- 
ture is determined by the balance between radiative heating and 
cooling processes: Because heating is mainly due to ionisation 
of hydrogen and helium, and because generally a significant con- 
tribution to cooling comes from line radiation of heavier ions, 
the models must become hotter with lower metallicities. During 
the course of evolution, the electron temperatures generally in- 
crease because the stellar photons become more energetic (see 
discussion in Sect. 3.1.2 above). At low metallicity the contribu- 
tion of radiation from collisionally excited metal ions is reduced, 
and line cooling is limited mainly to free-bound emission of hy- 
drogen and helium. Hence, cooling by expansion becomes more 
important for the energy balance, thereby limiting the electron 
temperature increase. This process is favoured by the larger ex- 
pansion rates of the low-metallicity models. See Sect. 3.1.4 for 
more details. 

The bottom panels of Figs. [5149] display both the dynamic 



and the equilibrium electron temperatures. The close correspon- 
dence between the dynamic and equilibrium electron tempera- 
tures seen in the metal-rich models during the horizontal evolu- 
tion across the HR diagram indicates that the gas heating time 
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Fig. 18. Mean [O m] electron tempera- 
tures, computed according to Eq. {T}, 
for all a = 3 nebular models as func- 
tion of the post-AGB age (left) and stel- 
lar effective temperature (right). In the 
right panel, the models are only plot- 
ted until maximum central star temper- 
ature (corresponding to a post-AGB age 
of -9800 years in the left panel) in 
order to avoid confusion due to over- 
lapping lines. The grey area depicts 
approximately the region occupied by 
the hydrodynamic models presented in 
|Paper I| and |Paper IVj which have all 
Z = Z GD , but various initial conditions 
and central stars. 
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the electron temperatures can be higher by up to 30 % in the 



metal-poorest case of Fig. 19 
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Fig. 19. Ratios between electron temperatures of equilibrium 
(eq) and dynamical (dyn) a = 3 models with different metal con- 
tents are plotted versus post-AGB age. The filled dot on each 
curve indicates the position of maximum stellar temperature. 
The sharp local minima of the temperature ratios occur during 
nebular recombination while the central star fades rapidly (cf. 
Figs.[3]and[8]l. See text for details. 



scales are much smaller than the stellar time scales for heating 



due to the ionisation of hydrogen and helium (see also Marten 
1995} . This statement holds independently of the metal content 
of the gas, and any deviations seen at later stages and/or for lower 
metallicities must be due to the dynamics. As expected, the tem- 
perature differences become largest for the models with the low- 
est metallicity, and hence lowest line cooling, and can reach up 
to about 10 000 K in extreme cases. 



Figure 19 gives a more detailed illustration of how the de- 
viations between dynamical and equilibrium models, measured 
by their mean 2+ electron temperatures, develop with time (or 
evolutionary stage of the central star). One sees that, in gen- 
eral, the differences increase with evolution, and they become 
large for the metal-poor models: Under equilibrium conditions 



During the fast luminosity drop of the central star between 
10000 and 11 000 years (cf. Figs. [3] and [8} one sees in Fig. [19 
that the mean temperature ratio eq/dyn becomes equal to or even 
slightly lower than unity. Line cooling increased by recombina- 
tion is not sufficient to explain this; rather we see a typical time- 
scale effect: The fading time scale of the star, -L/L, becomes for 
about 300 years as short as 400 years, which is quite compara- 
ble to the cooling time scale of the diluted nebular gas. Thus, it 
may well happen that the electron temperature in thermal equi- 
librium comes close to or even falls below non-equilibrium value 
because the latter cannot keep pace with the fading central star. 
After stellar fading is completed (at a;ll 000 years in Fig. 19 1, 
the previous differences of electron temperatures are restored 
quickly. 

A few comments concerning the halo, i.e. the still rather 
undisturbed but ionised AGB material ahead of the outer shock 
front, are in order here. Marten (1993a) [1995] ) showed that, be- 
cause of its rather low density, the halo gas is especially prone 
to non-equilibrium conditions, even with solar metal content. 
Generally one can say that the energy balance of the inner, 
denser halo regions of young PNe are controlled by line cooling, 
while in the outer, less dense regions expansion cooling prevails. 

Right after the passage of the ionisation front, the halo be- 
comes quite hot and cools then slowly down, with a (local) time 
scale ruled by the local density. Such temperature profiles can be 
seen for the metal-rich models in Figs. [6|j9] (lower leftmost pan- 
els) which are very compact and become optically thin quite late 
when the central star is already quite hot (T e ff 40000^-2 200 
K). 

The situation is different for the other models. They become 
optically thin earlier with cooler central stars (T e ff 33 000- 
35 000 K), and the halo is not so much heated in the first place. 
According to Marten ( 1995 1 the mean heating time scale of the 
halo gas can be comparable to or even larger than the stellar heat- 
ing time scale due to ionisation of hydrogen and helium. Under 
such condition the dynamical haloes remain always quite cool 
during the whole evolution, in contrast to the equilibrium con- 
dition. In general, the halo temperatures follow the trend of the 
main nebulae and become hotter with decreasing metal content. 
In the extreme cases the equilibrium electron temperatures ex- 
ceed locally 40 000 K (cf. Figs. Mf). 
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lonisation structure Although a PN might well be out of ther- 
mal equilibrium because of dynamics, ionisation equilibrium is 
generally rather well established. During the high-luminosity 
phase as the object crosses the HR diagram, ionisation dominates 
over recombination in the ionisation equations, and the nebu- 
lar ionisation quickly adjusts to the stellar ionising photon flux. 
Our hydrodynamical models are thus always close to ionisation 
equilibrium. An exception is only possible for a brief period dur- 
ing the rapid stellar fading when the ionisation time scale may 
become larger than the recombination time scale (e.g., of hy- 
drogen), and provided the latter is then comparable to or larger 
than the fading time scale of the star. Ionisation equilibrium is 
quickly re-installed later when the stellar luminosity evolution 
slows down on the white-dwarf cooling track. 



As an example we compare in Fig. 20 the ionisation proper 



ties of the dynamical and the thermally relaxed Zgd/10 model 
from Fig. [7] We selected this model because it is far evolved and 
has a very high degree of excitation because of the very hot and 
luminous central star. The metallicity chosen, Zgd/10, is typical 
for very metal-poor populations, e.g., in the Galactic halo and in 
some distant stellar systems. 

One sees in Fig. 20 that behind the shock (r 12 x 10 17 cm), 
i.e. in the nebula, the ionisation structures of both model types 
are very similar but not identical. The deviations seen in some 
cases are up to 10 to 15 % for main ionisation stages. These ob- 
vious deviations from a pure p/iofoionisation equilibrium are due 
to collisional ionisation whose contribution, albeit quite small, 
increases with electron temperature. Different electron temper- 
atures lead then consequentially also to small differences of the 
ionisation. 

The situation is demonstrated for hydrogen and helium (2nd 
and 3rd right panels in Fig. |20) . In the dynamical case, hydro- 
gen and helium are less ionised both in the halo and the nebula 
proper. The departures of the ionisation fraction ratios from unity 
follows closely the difference between, or ratio of, the electron 
temperatures. Only during the passage through the of the high- 
temperature shock region the ionisation of hydrogen and helium 
becomes temporarily larger than in equilibrium because of en- 
hanced collisional ionisation. 

Note also that for the particular, very diluted model shown 
in Fig. 20 where the central star is very hot and still quite lumi- 
nous, doubly-ionised oxygen which is so important for analysing 
nebular spectra is only a minority species. 

Because of its lower density, the halo, i.e, the matter ahead of 
the shock at r 14 x 10 17 cm, shows generally a larger degree 
of ionisation than the main nebula. For instance, the fraction of 
neutral hydrogen is about 10~ 4 in the halo, but about 10~ 3 in 
the nebula. Similar differences occur for the fraction of singly 
ionised helium: about 4 x 10" 3 in the halo and about 4 x 10 -2 in 
the nebula. The heavier elements show the same behaviour (cf. 
left panels of Fig. [20) for N, O, and Ne). 

It is also seen in Fig.|20](right panels) that the halo region is 
also prone to deviations from the ionisation/recombination equi- 
librium: the haloes of the dynamical models have lower degrees 
of ionisation than in the corresponding equilibrium cases (cf. 
nitrogen, oxygen, and neon in Fig.|20|i, indicating that the gas 
ionisation time scale exceeds the stellar ionisation time scale. 
The reason is the small fraction of high-energy photons which 
are able to sustain such a high ionisation. Thus, the ionisation 
time scale becomes the longer the higher the gas ionisation. 
Additionally, ionisation by electron collisions is more important 
in the cases of the very high halo temperatures. In equilibrium, 
the 5th ionisation (N, O, and Ne) appears to be the preferred 
one. After the halo gas becomes swallowed by the shock recom- 



bination brings the gas immediately back to a lower degree of 
ionisation very close to equilibrium conditions. 

The generally higher degree of ionisation within the haloes 
leads to a lower line cooling efficiency which in turn is respon- 
sible for the quite substantial temperature jumps across the neb- 
ula/halo boundaries, i.e. the shock, seen in the equilibrium mod- 
els (top right panel of Fig. 20 1. 



3.1.4. Line vs. expansion cooling 

Here we want to discuss in more detail the relevance of the dif- 
ferent cooling processes encountered in our models. The nebu- 
lar gas is heated by ionisation and cooled by line radiation from 
free-free, bound-free (recombination) transitions, and from col- 
lisionally excited atomic or ionic levels|^Additionally, the gas is 
subject to dynamical processes leading to cooling by expansion 
or heating by compression. 

In standard photoionisation modelling it is implicitly as- 
sumed that dynamics is unimportant, hence the thermal balance 
of the heated nebular gas is controlled by line radiation only. 
Since the most important coolants are metal ions, the cooling 
function decreases with metallicity until radiation from hydro- 
gen and helium prevails. Furthermore, radiation cooling depends 
on gas density squared, leading to a reduction of the cooling ef- 
ficiency also because the nebula expands with time. In contrast, 
cooling by expansion depends only on the gas velocity. 

Under adiabatic conditions, the total change of thermal en- 
ergy content of a volume element, U (in erg cm -3 ), with time is 
given by 



DU 



-p (V ■ v) - U (V ■ v), 



(2) 



where v is the flow velocity and p the thermal pressure (see Eq. 
All in |Marten & Szczerba|1997| l. The term p (V • v) is the usual 
local source (sink) term due to compression (expansion). The 
second term of the left-hand side of Eq. (|2]) accounts for the 
change of the volume of a gas parcel with time, i.e. a gas par- 
cel becomes expanded or compressed while streaming. 

It is more convenient to follow a mass element and rewrite 
Eq. $ using the thermal energy per mass, e = U/p (in erg g '), 
as 



P 



De 

d7 



-p(V-v). 



(3) 



Note that only in case of constant density, pDe/Dt - DU/Dt 
follows. 

Introducing p = nk^T e , with n being the total particle den- 
sity (ions + electrons), k B the Boltzmann constant, and T e the 
electron temperature, we get from Eq. Q in spherical coordi- 
nates 



De \dv 2v\ -. 

Ptt- =-nk B T e < — + — \ [erg cm s ]. 
Df \or r ' 



(4) 



Estimating possible dynamical effects demands knowledge 
of the density and the velocity structure of the object in question, 
both of which are difficult - if not impossible - to get observa- 
tionally. One may use the measured line width Vhwhm which, 
however, is an unsuitable choice: This velocity value is (i) much 
lower than the real expansion speed and (ii) is constant, (i.e. the 



3 Our time-dependent code treats heating by ionisation and cooling 
by line radiation separately, without considering net heating rates which 
hold only under the condition of strict ionisation equilibrium. 
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velocity gradient disappears), leading to a substantial underesti- 
mate of expansion cooling. The assumption of a constant elec- 
tron density may also introduce an uncertainty. It remains the 
likewise difficult choice of a - distance dependent - radius range. 
One could determine the radial extent of an object from a model, 
or select the outer edge of the nebula. In the latter case, an addi- 
tional systematic underestimate of the expansion cooling would 
result. 



The whole situation is illustrated in Fig. 21 for an evolved 
nebular model of the a = 3 sequence with Zgd/10 where the 
central star is at its maximum effective temperature (cf. Figs. [7] 
and 20 1. We selected this metallicity because we have previously 
seen that at this value hydrodynamical effects become significant 
for more evolved and diluted models. This figure allows both a 
detailed insight into the behaviour of the different cooling pro- 
cesses and an assessment of dynamical cooling. 

We employed the equilibrium temperature profile for esti- 
mating hydrodynamical effects since this is the case if one uses 
a photoionisation code. For this metallicity and ionisation stage 
line cooling is dominated by free-bound transitions throughout 
the entire nebula (Fig. 



21 



middle left). The largest contribu- 
tion to cooling from collisional excitation comes from neutral 
hydrogen, albeit its fraction is very small (« 10~ 3 ). Expansion 
contributes to the total cooling everywhere, except at the shock 
where the gas is heated considerably by shock compression (Fig. 



21 



middle right). Expansion cooling makes up for about 35% 
of the line cooling in the outer regions of the model immedi- 
ately behind the shock where the expansion is the fastest and the 
density the lowest. This extra cooling is responsible for the tem- 
perature difference of a few thousand degrees seen between the 
dyna mical and equilibrium model in the right top panel of Fig. 

m 



The bottom left panel of Fig. 21 compares the relative impor 



tance of hydrodynamical cooling for additional models of the 
a = 3 sequences. All these models are of the same age, which 
corresponds to maximum stellar temperature (cf. Fig. [7J. For 
these models, the contribution from hydrodynamics must be con- 
sidered if the metallicity drops below 1/3 to 1/10 of the Galactic 
disk (or solar) value. Furthermore, we also see that the hydro- 
dynamic cooling in the metal-rich models is restricted to outer 
regions. In their inner regions, i.e. the rims, a negative velocity 
gradient indicates compression which, at least partly, compen- 
sates the second term within the curly brackets of Eq. Q. This 
behaviour disappears with decreasing metallicity. One can also 
see that the shapes of the curves representing the ratios of to- 
tal cooling to line cooling correspond directly to the electron 
temperature differences between the dynamical and equilibrium 
models in Fig. [7] 



The bottom left panel of Fig. 21 also shows that for all these 



very advanced models the dynamics is important for the thermal 
balance of the halo regions: Expansion cooling is larger than line 
cooling for all shown models, although the former is quite small 
because of the low flow velocities and a slight compressional 
contribution from a tiny negative velocity gradient (see top left 
panel and right panel of the second row, both for r ^ 14 X 10 17 
cm). 



The right bottom panel of Fig. 21 illustrates the errors made 
when Eq. Q is used with various simplifications. Generally, the 
dynamical cooling is always underestimated, in one case locally 
up to a factor of ten! The best result is achieved by using the cor- 



4 It should be remarked that for the dynamical case all the cooling 
contributions seen in Fig.[2T|are slightly changed because of the differ- 
ent electron temperatures. 



rect density and radius variables of a model (provided they are 
available) together with the Vhwhm velocity. Since the electron 
temperature is fairly constant through the whole nebular struc- 
ture, an empirical mean value is a reasonable choice for all cases. 
The error figures found here are model-dependent and cannot be 
generalised. 

We conclude from the discussion in this section that any 
method to estimate empirically the contribution of dynamics to 
the total cooling function from an otherwise static model, as 
done by Stasihska et al. (2010 1 in their study of the metal-poor 
PN G135.9+55.9, will not provide a sound result. These authors 
used Eq. Q with the assumption p De/Dt = DU/Dt, but omitted 
the velocity gradient by setting u(r) = Vhwhm and r = R out (cf. 
right bottom panel of Fig. 21 1. - Only from full dynamical sim- 



ulations one can get hold of the correct contribution of dynamics 
to the heating/cooling balance of a nebular model. 

3.1.5. Line strengths 

The temperature differences between dynamical and equilibrium 
models have, of course, consequences for the line emission, es- 
pecially for that of the collisionally excited lines, as has been 
already demonstrated by Schonberner et al. (2005c). The effect 
is again illustrated in Fig. 22 for three lines of oxygen, in the UV 



(O iv] A 1402 A), the optical ([O m] A 5007 A), and the IR region 
([Oiv] A26 /im). The trend is as expected: 

- In equilibrium the models show a tendency to have always 
higher line strengths, in line with the temperature difference 
to the dynamical case. 

- The differences increase with evolution (or stellar tempera- 
ture) and decreasing metallicity; they can, however, be to- 
tally neglected for more normal or metal-rich compositions 
(cf. Fig. 22 left bottom panel). 

- The relative differences of line strengths between both types 
of models are the largest for UV lines and the smallest in 
the infrared. The optical lines behave intermediately. At the 
lowest metallicities considered here (Zgd/100), Oiv] A 1402 
A can be off by about 250 % if thermal equilibrium is as- 
sumed. The corresponding value for [O iv] A 26 /im is only 
about 30 % (cf. Fig. [22] right bottom panel). 

Note that the main reason for the differences of line strengths 
between equilibrium and dynamical models is the electron tem- 
perature because we found only very small deviations from ion- 
isation equilibrium (cf. Sect. 3.1.3 above). We repeat that the 



discrepancies found here are highly dependent on structure and 
expansion properties of the model and should be regarded as in- 
dicative only. 



We see also in Fig. 22 that the expected changes of line 
strengths with the corresponding abundances is partly compen- 
sated by the electron temperatures. For instance, at medium stel- 
lar temperatures when 2+ is the main ionisation stage and the 
electron temperatures increase with decreasing Z (cf. Fig. 18 i, 
the line strength of 5007 A decreases roughly with 0.5e(O) only. 
At high excitations (T e g ^ 100000 K) when most oxygen is 
3+ , the electron temperatures of very metal -poor (dynamical) 
models become nearly independent of metallicity (Fig. 18 1, and 
hence the strength of, e.g., 26 fim varies proportional to e(O). 

Our metal-poor models are well suited to address the ques- 
tion of how one can determine the chemical composition of 
metal-poor PNe if conventional plasma diagnostics is not fea- 
sible because of lack of suitable lines, as is the case for 
PN G135.9+55.9. One is then left with photoionisation mod- 
elling for constraining, e.g., the oxygen abundance, as has been 
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Fig. 22. Top panels: run of line strengths (relative to HB = 100) of three collisionally excited oxygen lines, Oiv] A 1402 A, 
[Om] A 5007 A, and [Oiv] A26 pm, with stellar temperature for three a — 3 model sequences, viz. Z = Zgd (left), Z - Zgd/10 
(middle), and Z = Zgd/100 (right). The evolution is only shown until maximum T e ff in order to avoid confusion. Note the different 
ordinate ranges. Bottom panels: Line strength ratios L eq /L^ yn . The small scatter seen for the A 1402 A ratio in the metal-rich model 
is due to numerics. 



shown extensively and illuminatingly by Pequignot & Tsamis 
(2005). According to these authors, the oxygen abundance of 
PN G 135. 9+55. 9 is likely not below 1/30 of the solar value, as- 
suming a stellar temperature of a 130 000 K^] 

However, for a highly-excited PN like PN G135.9+55.9 the 
oxygen content can best be estimated by using, e.g., the 26 pm 
line whose strength depends only weakly on the electron temper- 
ature and its uncertainties (cf. Fig. 22 1. The optical oxygen line 
at 5007 A cannot be recommended because its higher sensitivity 
to the electron temperature, especially if the stellar temperature 
is above 100000 K. Ultraviolet lines of highly ionised elements, 
such as O iv] A 1402 A, should not be used at all for abundance 
studies because these lines depend strongly on electron temper- 
atures at nearly all stellar temperatures. 

We conclude from this discussion that the use of photoioni- 
sation models assuming thermal equilibrium is not adequate for 
objects with a diluted nebula and a very hot central star (i.e. for 
highly excited objects as represented by our models) and with 
reduced metallicity because of incorrect physical assumptions. 
Note that departures from thermal equilibrium depend only on 
the gas properties regardless of the nature of the central ionising 
object. Of course, this will become evident only if the restriction 
to one wavelength region (e.g. the optical) is relaxed. Our mod- 
els show that, for UV lines, departures from equilibrium may be- 
come a problem already for Z ^ Zgd/3. They suggest also that 
the neglect, or a too simple consideration, of dynamical effects in 
the modelling of a metal-poor and highly excited PN, as done in 
the recent study of PN G135.9+55.9 by |Stasiriska et al.| ( [20T0] l, 
becomes problematic. 

Based partly on new spectroscopic material, Sandin et al 



( |2010| l performed an abundance study of PN G135.9+55.9 tak- 
ing full advantage of the dynamical models discussed here. For 
the details, the reader is referred to the cited paper. 

3.1.6. The Milky Way halo PNe 

In this context it is interesting to see how the known Galactic 
halo PNe fit into our grids of mean electron temperatures vs. stel- 



lar effective temperatures or ages of the simple hydrodynamical 
a — 3 models as shown in Fig. 18 and whether hydrodynamical 
effects are discernable. As ages of the halo PNe are not known, 
and their stellar temperatures are too uncertain, we used the neb- 
ular excitation as given by the relative strength of He n A 4686 
A for a proxy of the evolutionary status of both the objects and 
models. 



The result is seen in Fig. 23 where the mean [O in] tem- 
peratures are now plotted against L(4686)/L(H/3) instead. We 
compared our models (also in equilibrium) with both the elec- 
tron temperatures from detailed nebular analyses by means of 
photoionisation models as provided by Howard et al.| ( 1997} 
and with the corresponding temperatures derived directly from 
the observed temperature sensitive [Om] line ratios 648 and 
DdDm-1 are omitted because they do not show a He n A 4686 A 
line. 

Neglecting for the moment M2-29, the following emerges 
from Fig. 23 



Compared with the simulations, the observed mean electron 
temperatures suggest nebular abundances between Zgd/3 < 
Z < Zgd/30, which is in reasonable agreement with the spec- 
troscopic analyses in the literature (cf. Table 4 in Howard et 
|al.|1997| ). 

On average, discrepancies between observed and modelled 
temperatures increase with excitation: they are less than 
1000 K at L(4686)/L(H/3) a 0.1, but up to 3500 K at the 
highest excitation, L(4686)/L(HyS) - 1 PI This is consistent 
with the predictions of our model sequences if one neglects 
the (unexplainable) fact t hat for the 3 objects with medium 
4686 A line strengths the Howard et al.| ( |1997| temperatures 
are below the temperatures derived from the observed line 



5 The analysis of Pequignot & Tsamis 1 2005 I is entirely based on 
optical lines ! 



6 The electron temperatures which we derived from the observed 
[O m] line ratios listed in Table 2 of Howard et al. ( 1997 ) differ in some 
cases, for unknown reasons, from the temperatures given in the same 
table. 

7 The trend seen for the three highly excited objects is most likely 
not real. Firstly, the relative importance of expansion cooling depends 
on the object's density and velocity structure, and secondly, the electron 
temperature based on a photoionisation model is subject to the imposed 
fit criteria. 
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strenghts, while they are above for the 3 objects with the 
highest excitations. 

Because abundances derived from collisionally excited lines de- 
pend strongly on the electron temperature, one is tempted to say 
that the use of standard photoionisation models cannot be rec- 
ommended for cases where non-equilibrium (i.e. dynamical) ef- 
fects become important. This statement may become applicable 
already at metallicities below about one third solar, for very di- 
luted and rather fast expanding nebulae. 

For M2-29 with its only moderate excitation the discrepancy 
between observed (from Rom) and modelled temperatures is the 
largest, ^8000 K, and at disturbing variance with the model pre- 
dictions. First of all, the photoionisation models of |Howard et| 



al. 



are based on the incorrect [Om] /14363 A line strength of 
15 (HjS = 100) which is the undereddened value (cf. |Pena &| 
Torres-Peimbert| 1 99 1 1 Table 2A). We used in Fig. 23 the dered- 
dened value of 19.5 instead, leading to an electron tempera- 
ture of 24 000 K which is then the highest in the |Howard et al.| 
1997 ) sample of Galact ic halo PNe. The photoionisation model 



of |Howard et al.| ( |1997] > predicts a /14363 A line strength of 11 
only instead of the assumed value 15, which would correspond 
to an electron temperature of about 20 000 K. 

Because of the high electron temperature, the plasma diag- 
nostics of Pena & Torres-Peimbert ( 1991 ) leads to a lower metal - 



licity for M2-29 than claimed by |Howard et aL ( 1997] >: a metal 
depletion by factors between 20 and 30. Still, the problem of 
M2-29's too high electron temperature remains a mystery: At its 
moderate excitation this high electron temperature would sug- 
gests a metallicity well below Zgd/100, which is not observed. 
In this context it is interesting to note that Haj duk et al.| ( 2008 ) 
claimed recently that the nucleus of M2-29 consists of a binary 
system with an eclipsing disk. 

3.2. The a - 2 sequences 



Recently we showed in Pap er IV| that the bright cutoff of the 
PNLF can be explained by nearly optically thick nebulae around 
central stars with masses slightly above 0.6 M©. However, in 
Paper IV we only assumed a metal content that is typical for 
Galactic disk objects. We thus continue our work with inves- 
tigating the properties of model nebulae around more mas- 
sive (and much faster evolving) central stars, 0.696 M© and 
0.625 M©, and with different metal contents, 3Zgd, Zgd, Zgd/3, 
and Zgd/10. This is also the mass range of interest in observa- 
tions of very distant stellar populations because there only the 
most luminous nebulae are accessible for more detailed stud- 
ies (see, e.g., [Arnabold i et al.|[2008j ). All initial models have 
the same constant-outflow conditions, viz. p oc r~ 2 , M ag b = 10~ 4 
M©yr _I , and o ag b = 15 kms" 1 (cf. Sect. El. The metallicity de- 
pendent central-star wind model used is tne same as introduced 
in Sect. [2] For these always rather compact and dense models 
the dynamic contribution to the cooling function remains always 
negligible, also at lower metallicities. 



3.2.1. Expansion properties 

We start with a discussion of the expansion properties of models 

con- 



around a 0.625 M© central-star. For this purpose, Fig. 24 



tains four snapshots, each taken at the moment of maximum 
[O in] 5007 A luminosity of the respective sequence, showing the 
model structures and the corresponding profiles of the 5007 A 
line, both spatially resolved and volume integrated. The four se- 



quences started from exactly the same initial configuration (see 
above) but with different metallicities: 3Zgd, Zgd, Zgd/3, and 
Zgd/10. Relevant model parameters are listed in Table[3] 

The model structures of all sequences shown are very simi- 
lar, with a pronounced shell/rim morphology. The mass fractions 
contained in the rims are, however, quite different (Table[3j Cols. 
1 1 and 12) and reflect the dependences of stellar wind power and 
shock speed on metallicity: 

1 . For a given time, the total nebular mass, M; on (/? out ), increases 
strongly with decreasing metal content because of the in- 
creasing shock speed, R out (Col. 7 of Table^. 

2. The case for the rim is more complicated: In general, the de- 
creasing wind power leads also to a decreasing bubble pres- 
sure, and hence also to smaller rim masses. This trend is in- 
verted for the metal-poorest sequence because here the bub- 
ble becomes so small that its pressure increases again some- 
what, despite of the very weak stellar wind. 

3. In all cases, however, the nebular mass fraction contained 
in the rim remains rather small as long as the shell is not 
recombining (see below)! 

A remark concerning the terminology "young" planetary 
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nebula is in order here. The bright models shown in Fig 
and whose parameters are listed in Tableware all very compact 
(^out < 0.07 pc) and less than 2500 years old. Objects with such 
properties are usually considered to be young PNe. However, 
because of the higher mass (0.625 M©) of the central stars these 
systems are already quite evolved with hot stars of effective tem- 
peratures ranging between 55 000 and 120 000 Kelvin! The pres- 
ence of a rim is the result of this advanced stage of evolution. 
Thus, in terms of evolution such objects must be considered at 
least as being middle-aged, but not "young". The Milky Way 
object NGC 7027 is an even more extreme example of this kind: 
Because its central star is already close to its maximum temper- 
ature, the term "old" is a better designation for this PN (for more 
details on this object, see |Schonberner e"t~a l. 2005b Paper III 
hereafter). 

Concerning their general expansion properties the models 
with a = 2 behave (qualitatively) in the same way as the a — 3 
models discussed extensively in Sect. 3.1 A lower metallicity 



provokes higher electron temperatures and faster expansion, but 
a lower pressure exerted by the shocked stellar wind (top panels 
of Fig. 24 1. Consequently, the velocity difference between the 
outer and the inner nebular boundaries increases with decreas- 
ing metallicity as well. This effect is clearly seen in the middle 
panels of Fig. 24 where the computed, spatially resolved [O m] 



5007 A line profiles of the models are shown: At lower Z, the 
fast shell and the slow rim matter are spectroscopically resolved, 
and the entire profile becomes wider. In the Z = Zgd case both 
velocities are already too similar and cannot be separated. Note 
that in all cases shown (except for Z = 3Zgd where the rim mat- 
ter expands faster than the shell matter) the extreme values of the 
profile derivative, measured at the outer flanks of the profile, un- 
derestimates only slightly the post-shock velocity, V" post , i.e. the 
gas velocity immediately behind R out (see also the discussion in 
ICorradi et aLpOOT) . 

Spatially resolved line profiles are not available for PNe in 
extragalactic systems, and we thus show the corresponding spa- 



tially integrated profiles in the bottom panels of Fig. 24 as well 



Despite of the different velocity and density profiles of the mod- 
els caused by their different metallicities and which result in 
quite different spatially resolved line profiles (as is discussed 
above and shown in the middle panels of Fig. |24|), all four spa- 
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Fig. 24. Composite figure showing one snapshot taken at maximum 5007 A brightness along each of the four 0.625 M0 sequences 
(a = 2), with the metallicities indicated over the panels. Details of the models are listed in Table [3] Top panels: structures of the 
models: run of ion density (thick solid line), flow velocity (thin solid line), and electron temperature (dotted line) with radius. 
Densities and temperatures are normalised and must be multiplied with the factors given in each panel. The radial positions, R out , of 
the nebular outer shock and of the rim, R r [ m are also marked in each panel. The central star is at r — cm. The electron temperature 
spikes seen at larger radii are caused by the ionisation fronts, and the jumps at the inner nebular edges mark the positions of the 
contact discontinuity/conduction fronts. Middle panels: normalised spatially resolved [Om] 5007 A and thermally broadened line 
profiles (dotted line) as seen through a central numerical aperture with a diameter of 1 x 10 16 cm (corresponding to 0'.'67 at 1 kpc), 
and additionally broadened by a Gaussian of 10 kms 1 FWHM (solid line) in order to simulate a typical instrumental broadening. 
The maximum velocity behind the outer shock at R out , Vp OSt , is indicated by vertical dashed lines. Bottom panels: [O in] 5007 A line 
profiles as integrated over the whole nebula (dotted line), and additionally broadened by a Gaussian of 10 kms~' FWHM (solid 
line). Two vertical dashed lines in each panel indicate the HWHM and the shock velocity, R out , respectively. 



tially integrated profiles are very similar, but with much smaller 
line widths than found for their spatially resolved counterparts. 

The expansion rates as they follow from the line widths, 
Vhwhm, vary only between 18 and 20 kms~'and are well below 
y post (-24 ... 30 kms" 1 ), and only about half of R om (^33 . . .40 
km s~') - the true expansion rates (cf. Table[3] Cols. 7, 8, and 9)! 

The surprisingly low Vhwhm can be understood as follows: 
By integrating over the nebular volume the fraction of nebu- 
lar gas with the highest line-of-sight velocities is small and is 
only responsible for the faint wings of the total line profile. The 
high-speed matter immediately behind the outer shock is only 
detectable either by means of spatially resolved high-resolution 
spectroscopy ( |Corradi et al.|2007 1, which is not possible for dis- 
tant objects, or by measuring the width of the integrated profile 



close to the bottom of the profile. This, however demands a suf- 
ficiently high signal-to-noise ratio. 

The case is illustrated in more detail in Fig. [25] where the 
time evolution of various mass fractions of the model nebula 
with Z = Zgd around the 0.625 M0 central star is provided. One 
sees a large temporal variation of the total ionised nebular mass: 
a mass increase due to ionisation and the expanding shock which 
swallows the upstream AGB matter and adds it to the nebula, 
followed then by a rapid decline when recombination destroys 
nearly the entire shell after about 3500 years. For the density 
and velocity profiles of this model at an age of 1694 years, see 
Fig. [24] (second top panel). 



Figure 25 can be considered as a tool for interpreting the 
kinematics of a PN from spatially-integrated line profiles: Thick 
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Table 3. Relevant parameters of the four a — 2 models with central-star masses 0.625 Mq, as depicted in Fig. 24 Columns 2-A 
describe the central star (post-AGB age, luminosity, and effective temperature), and Cols. 5-12 the nebular properties (maximum 
[O in] -brightness, rim/shock positions, shock propagation, post-shock velocity, velocities at half-width-half-maximum (HWHM) 
and half-width- 10% -maximum (HW10%M) of the spatially integrated [O ra] 5007 A line profile, the mean 2+ electron temperature, 
the mass of the rim, and the total ionised nebular mass). The initial parameters of the circumstellar envelopes are M ag b = 10~ 4 
Mq yr -1 and i> ag b = 15 kms~' in all cases. 



Metall. Z 

(1) 


t 

[yr] 

(2) 


L 
[L©] 
(3) 


T e e 
LKJ 
(4) 


M(5007) 
[mag] 
(5) 


^rim/^out 

[10 17 cm] 
(6) 


[kms -1 ] 
(7) 


Impost 

[kms- 1 ] 
(8) 


,,[0111] ,,,[0111] 
'HWHM' V HW10%M 

[kms- 1 ] 
(9) 


[K] 
(10) 


Mon(#rim) a 

[Mq] 
(11) 


Mon(^ out) a 
[M©] 
(12) 


3Zqd 


2392 


6643 


1 17 643 


-4.66 


1.5/2.1 


33 


24.1 


20.3/ 26.3 b 


8 070 


0.043 (0.017) 


0.174 (0.109) 


Zgd 


1694 


7507 


73 338 


-4.25 


0.8/1.5 


35 


25.7 


18.6/25.1 


9 553 


0.012 


0.131 


Zgd/3 


1460 


7667 


60917 


-3.80 


0.6/1.3 


37 


27.6 


18.5/26.7 


12 080 


0.007 (0.011) 


0.119(0.153) 


Zgd/10 


1324 


7736 


54 296 


-3.08 


0.4/1.2 


40 


29.6 


18.0/27.9 


15 198 


0.012 (0.013) 


0.115 (0.169) 



a Mass values in parentheses refer to the age of the Zqd model, / = 1694 yr. 



b Vhwio%m > Vp OSt because at this phase of advanced evolution the rim matter expands faster than the shell (cf. top left panel of Fig. 24 1 



lines (solid for [Ora] 5007 A and dashed for H/3) indicate the 
mass fractions where the real gas velocities are lower than the 
value deduced from the corresponding half width of the total 
(volume integrated) line profile. This matter contributes fully to 
the line emission within the halfwidth value, next to matter with 
projected velocities < Vhwhm- The important velocity informa- 
tion from the faster moving gas, the mass fraction of which may 
even exceed half of the total ionised mass, can only be retrieved 
from the outer wings of the line profile. 

Obviously the only method to get hold of the signature of 
the fastest expanding nebular parts is to measure the line width 



at a very low (ideally zero) intensity, as proposed by Dopita et 
[aT] ( 1985"| |1988| l who suggested to measure line widths at the 
10% level from maximum (see also discussion in Sect. |4.4| i. One 
sees indeed from Fig. 25 that such a choice, Vhwio%m, ensures 



that the mass contained in the outermost shell in which the (real) 
gas velocities are higher than Vhwio%m becomes very small, i.e. 
Vhwio%m is close to the post-shock velocity y post (TableBl Cols. 
8 and 9). 

The situation becomes easier for this particular model when 
recombination sets in after 3400 years of evolution: The shell 
disappears nearly completely, and because the velocity within 
the rim is roughly constant, Vhwhm represents, in the particular 
case shown here, the overall PN kinematics well (Fig. 25 1. The 
now rather low PN mass (<0.1 M©) corresponds roughly to that 
of the rim formed by the stellar wind action during the previous 
high-luminosity part of evolution. 

The complete velocity evolution of three 0.625 M© se- 
quences is rendered in Fig. 26 Plotted are three velocities, R out , 
Vpost, and Vhwhm as a function of the post-AGB age. We see that 
in general the difference between the measurable Vhwhm and the 
real expansion speed, R out , is quite large: about a factor of two 
during the [Om] bright phase of evolution (e.g., filled circles in 
the figure). During the early, optically thick stage, the velocities 
increase rapidly with time in all cases. The Vhwhm velocities 
start at very low values, well below the AGB wind velocity of 15 
kms" 1 , because early in the PN evolution when the central star 
is not very hot, 2+ is restricted to the inner, slowly expanding 
matter in or around the rim only (cf. Fig. 24 top panels). 

As we have already seen for the a = 3 sequences, Vhwhm 
is quite insensitive to the metal content for the whole evolution, 



in comparison to the shock speed R mt (cf. Fig. 13 1. However, 
the post-shock gas velocity, V post reflects the metal dependence 
of the shock speed quite well (Fig. [26), but it requires accurate 



measurements of the faint wings of the integrated line profile, as 
discussed above. 

The velocity bumps (Vp OSt and R oat only) seen close to 
t ot 1000 yr indicate that at this time the models become optically 
thin for Lyman continuum photons and enter the 'champagne' 
phase of expansion during which the shock speed depends only 
on electron temperature (sound speed) and radial density gradi- 
ent (see [Paper II) . This transition is not seen in Vhwhm because 
the outer nebular region behind the shock contains only little 2+ 
at this time. Later, at t =s 1200 yr, also the outer nebular regions 
start to contribute significantly to [O m] line profiles, and the ve- 
locity increase flattens out. The deceleration of the shock speed, 
R out , at higher ages, t > 3000 yr, is due to recombination behind 
the shock, leading there to much lower electron temperatures and 
hence lower pressure. 

Of course, the maximum amount of [Om] emission de- 
creases with the oxygen content (see Table |3). This table also 
shows that this maximum emission value occurs the earlier (i.e. 
at lower stellar temperatures) the lower the oxygen abundance 
is (cf. filled circles in Fig. 26 1. The reason for this behaviour 
is twofold: (i) the + ionising continuum becomes optically thin 
earlier for lower abundances, and (ii) the number of photons ion- 
ising + (A < 353 A) decreases rapidly with T e ff 

We note also that the expansion speeds do not reach the very 
high levels of the models shown in Figs. 12 and 13 It has noth- 



ing to do with the different central-star masses (0.625 vs. 0.595 
M©). Rather it is caused by the different density profiles of the 
circumstellar matter. We recall that the nebular models shown in 
the Figs. 12 and 13 have initial models with power-law density 
with a — 3, while the sequences discussed 
2, which is the signature of an AGB 



profiles, p{r) cc r 
in this section have a 
wind with fixed speed and mass-loss rate. For a — 2, the shock 
speed relative to the upstream flow is considerably lower than 
for the a — 3 case (see Fig. 7 in Paper II) . The modest shock 
acceleration during the (optically thin) 'champagne' phase seen 



in Fig. 26 is due (i) to the rise of the electron temperature as the 
star becomes hotter, and (ii) to the weak acceleration of the AGB 
wind (= upstream flow) after the passage of the ionisation front. 



8 The stellar radiation field is assumed to be independent of the cho- 
sen metallicity. 
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3.2.2. The [Om] brightness 

The [O m] emissiv i ty sho uld depend directly on the oxygen con- 



tent itself. |Jacoby 



1989) hat the 5007 A lu- 

minosity changes with roughly the square root of the abundance 
only. This more modest dependence is explained by two effects 
which compensate each other partly: The lower oxygen abun- 
dance leads to a higher electron temperature, which in turn en- 
hances the emission of collisionally excited lines. Based on op- 
tically thick photoionisation models, Dopita et al. ( 1992| l found 
that there exists a maximum efficiency of converting stellar ra- 
diation into [O m] 5007 A line emission at about twice the solar 
oxygen abundance]^] 

Since our a - 2 sequences (0.696 M0 and 0.625 M0) cover 
the range about the bright 5007 A cut-off of the PNLF, we used 
them also to check the variation of the 5007 A luminosity with 
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how the 



metallicity Z^For this purpose we show in Fig 
maximum 5007 A luminosity, relative to the bolometric lumi- 
nosity of the central star (i.e. the efficiency of converting stellar 
UV radiation into [O m] 5007 A line emission), varies with Z, 
and compare the outcome of our hydrodynamical models with 
predictions of (optically thick) photoionisation models of Dopita 
et al. ( 1992 1. We used their Eq. (4.1) which is a fit to their model 
predictions, corrected for the fact that Dopi ta et al.| used a higher 
solar oxygen abundance (8.84 on the usual scale). 

The maximum efficiency of converting stellar radiation into 
[0 111] 5007 A line emission is 0. 1 1 for the (optically thick) 0.696 
M0 models with the Galactic disk composition used here (i.e. 
e(O) = 8.65). The corresponding stellar effective temperature 
is about 110000 K. The maximum efficiency decreases with 
Z and occurs at progressively lower stellar temperatures, i.e. 
r eff a 136000 K for Z = 3Z GD , T eS 101 000 K for Z GD /3 
and ^95 000 K for Zgd/10. Our models confirm that this ef- 
ficiency decrease is weaker than expected because the reduced 
oxygen abundance is partly compensated for by an increase 
of the mean electron temperature: from (T e ) = 8170K (6Zgd) 
over 9500 K (3Z GD ), 1 1 440 K (Z GD ), and 13 980 K (Z GD /3) to 
16650K (Z G d/10). According to our models, the general Z 05 
dependence of the 5007 A conversion efficiency quoted in the 
literature is inaccurate: it holds only for PNe with a metal con- 
tent around Z GD /3 (see Fig.^] 17 ] 

Our models are, however, partly discrepant to the results of 
Dopit a et al.| ( [T9 92): For metal contents below the solar one the 
agreement between our hydrodynamic models and the Dopita et 
IT] fit is reasonably good, but according to their Eq. (4.1) the 
5007 A conversion efficiency reaches a maximum at Z — 3 Z G q 
not found with our models (cf. Fig. |27j ). Our models are rather 
suggesting a levelling off at extremely high oxygen abundances 
instead. We do not know the reason for this discrepancy, since 
our hydrodynamic models, being still quite dense and optically 
thick, are also in ionisation and thermal equilibrium. However, 
our hydrodynamical sequences mimic consistently the evolution 
of bright nebulae systems as realistic as possible, whereas the pa- 
rameter space of the Dopita et al. models is limited and does not 
cover the high stellar luminosities necessary for PNe close to the 



9 This statement refers to the old solar oxygen abundance of e(0)o = 
8.84. 

10 Note that Z refers to an abundance pattern scaled to those of the 
Galactic disk PNe. Thus, we always have Z oxy5C11 oc Z. 

1 1 We recall that the metallicity discussion refers only to the circum- 
stellar matter. The central-star models (and the radiation fields) used are 
the same in all four sequences and are based on evolutionary computa- 
tions with about solar initial composition (Blocker 1995} . 



luminosity function cut-off. Therefore, we caution that |Dopita et| 
aT]s fit formula for estimating the metallicity dependence of the 
[O m] cut-off brightness might lead to erroneous conclusions if 
it is extrapolated into the metal-rich domain (cf., e.g., |Ciar dullo 
et al.|2002| Fig. 5 therein). 

The models around the 0.625 M0 star become optically thin 
and therefore have always lower conversion efficiencies. A more 
detailed description how the [Om] brightness depends on model 
parameters and/or evolutionary stages can be found in |Paper fV| 
although only for the Galactic disk chemical composition. 



4. PNe in distant stellar systems 



We have already shown at different occasions, viz. |Paper II 



Paper III] |Steffen & S chonberner (2006), and recently in 



Paper IV that the basic physical model used in our simulation 



of the planetary nebula evolution is very successful in explain- 
ing basic observed structures of PNe. A further important test of 
the quality of our hydrodynamical models is to answer the ques- 
tion if and how well the models are also able to explain global 
properties of PNe samples in stellar systems so distant that a PN 
cannot be spatially resolved. Very useful for testing are basic dia- 
grams as those involving characteristic line ratios and expansion 
velocities measured from the Doppler width of bright emission 
lines. 



4. 1 . Line-ratio diagrams 

We begin with a diagram which combines the line ratio 
L(5007)/L(Ha + [Nn]) with the brightness in [Om], M(5007). 
As an example, Fig. [28] compares the data of two Local Group 
gal axies, M31 ( |Ciardullo et al.||2002| i and M33 ( |Magrini et 
|aLl|2"000), with our model predictions. The observations oc- 
cupy a wedge-like region in this diagram and are, according to 



Herrmann et al. ( 2008 1, roughly bounded from below by a line 

log{L(5007)/L(Ha; + [Nii])} = -0.37M(5007) - 1.16, 

and from above by a horizontal line at 

L(5007)/L(Ha+[Nn]) ^4. 

The lower boundary is nearly perfectly matched by all of 
our models during their early, optically-thick evolution towards 
maximum 5007 A emission, virtually independently of their pa- 
rameters and chemical composition. This boundary is not very 
well populated because this optically-thick phase of evolution is 



rather short (cf. Paper IV) Fig. 15 therein). The maximum reach- 



able brightness in M(5007) depends, of course, on the central 
star mass, the metallicity, and the thick/thin transition. Models 
which remain optically thick (i.e. the ones with the 0.696 Mq 
central stars) evolve then backwards a little bit above this lower 
boundary. Models which become optically thin turn then further 
upwards until L(5007)/L(Ha + [Nn]) ^ 4 is reached and stay 
there until the central star fades: The line ratio decrease after 
the 'kink' is the result of the rapid luminosity drop of the central 
star as hydrogen burning ceases, with a corresponding increase 
of the [Nn] lines because of a temporary recombination to N + . 
A slow increase of this line ratio occurs only at fainter mag- 
nitudes (M(5007) ^ 0) because of re-ionisation during the very 
slow early white-dwarf phase. Because of the fast model evo- 
lution between the 'kink' and the re-ionisation stage, it is very 
unlikely that observed objects are in this particular phase. 

We can therefore state that the wedge-shaped region of Fig. 
28 is populated by PNe which are evolving along the horizontal 
part of the HR diagram through maximum [O m] brightness until 
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the central star begins to fade. Nebulae around the most massive 
central stars, ^ 0.63 Mq, stay always rather close to the lower 
boundary because they remain optically thick or become only 
modestly optically thin. Optically thin (or partly thin) models 
around central stars of ^0.63 Mq to well below 0.59 Mq which 
are on their evolution across the HR diagram fill the entire region 
between the lower and upper boundary and towards the right 
until faint magnitudes. 

A similar diagram can be made for 5007/(Ha + [N n) vs. 
M(Ha + [N it]) and is shown in Fig. 29 In this diagram the PNe 
remain bright in (Ha + [N n]) during the fading of the central 
stars towards white-dwarf configurations because of recombina- 
tion and a temporary dominance of the [Nn] lines: The models 
are now in Ha + [N n] brighter than in [O in] 5007 A by up to 1 .5 
mag. Hence, our hydrodynamical models suggest that the deficit 
of objects in the lower right corner of this diagram is not due to 
sample incompleteness (as proposed by Ciardullo 2010) but is 
rather a consequence of central star evolution: The early white- 
dwarf phase where the nebula re-ionises limits the observability 
of PNe in distant systems if observed in Ha + [N n]. 

Another useful diagram for a comparison between data and 
the predictions of our simulations is that presented in Richer 
(2006 Fig. 1 therein). It combines the nebular ionsation struc 
ture, as measured by [O m] A50Q1 A and [O n] A3727 A, with 
the strength of Hen /14686 A, which is essentially an indicator 
of the stellar surface (effective) temperature. An adaptation of 
Richer s Fig. 1 is shown in Fig. 30 with the predictions of our 



evolutionary models included. This diagram is nearly indepen- 
dent of metallicity, as is demonstrated by the three 0.625 Mq 
sequences with different Z, but contains only objects with rather 
hot and thus evolved central stars because /14686 A must be mea- 
surable. 

The interpretation of this diagram, based on our evolutionary 
models, is as follows: 



- The right border, at L(4686)/L(H/3) ^ 1.0, is occupied by 
highly ionised objects for which the emitting volumes of Hy8 
and He n /14686 A are about equal. When recombination oc- 
curs, this ratio becomes smaller, as is seen in our models. 
At the same time, the strength of the [O n] doublet increases 
considerably, more than the /15007 A line, because the mod- 
els become optically thick. These objects close to the right 
border have very hot and still luminous central stars, produc- 
ing high nebular excitation, but at the same time also strong 
[O n] emission, and [N n] emission as well. 

- The upper limit L(5007)/[L(3726) + L(3729)] 100 is ex- 
plained either from above by optical thick sequences with 
massive central stars («0.65 Mq), or from below by se- 
quences with less massive central stars whose models be- 
come optically thin during the course of evolution. Objects 
with massive central stars with S 0.65 Mq are obviously not 
existent (cf. also Fig. 28 i. 

- A lower boundary for L(5007)/[L(3726) + L(3729)] is pro- 
vided by optically thick models around central stars of about 
0.6 Mq (the three outliers neglected). 



The whole observed sample can be interpreted to consist of 
objects with rather highly ionised/excited optically thick or thin 
nebular envelopes, as dependent on their evolutionary stage, with 
central stars with masses between about 0.6 and 0.65 Mq. 



4.2. Expansion properties 

In distant stellar systems one cannot spatially resolve the pro- 
files of emission lines, and information on expansion properties 
can only be gained from measuring the halfwidth of the pro- 
file. Although relevant studies are rare, the following picture 
1985 1988) showed for a large sam- 



emerged: Dopita et al. 



pie of Magellanic Cloud objects that there exists an evolution 
of nebular expansion rates, in terms of line widths, from low 
to higher values with nebular excitation class or stellar effec- 
tive temperature. Richer (2006 2007) reported, based on an ex- 
tensive spectroscopic study of PNe in the Local Group galax- 
ies, that there is virtually no dependence of the mean expansion 
properties, as measured from the halfwidth of the 5007 A [O m] 
line, Vhwhm(5007), on the metallicity of the parent stellar pop- 
ulation. This halfwidth varies between about 8 and 30 kms -1 , 
with a clear accumulation around 18 km s -1 for the brightest ob- 
jects (|Richer|2006| |2007[ ). For the Galactic bulge PNe, |Gesicki| 
|& Zijlstra||2005| found also a mean V H whm(5007) of 18 kmFT 
This obvious observed invariance of PNe emission line 
widths with metallicity seems to contradict strongly the predic- 
tion of our hydrodynamical models that expansion speeds gen- 
erally increase with decreasing metal content. We thus used the 
sequences discussed in the previous sections for clarifying the 
situation and interpreting existing observations. 

We start with showing the evolution of our models in the 
Vhwhm-^(5007) plane in Fig. [31] which is adapted from Fig. 4 
in |Richer| ( |2006| ) and contains also the same data from dif- 
ferent galaxies and/or parts of them. All the relevant data are 
now also available from |Richer e t al. (2010b). The theoretical 
tracks of our two most massive central stars (0.696 and 0.625 
Mq) surrounded with PNe models with different metallicities 
(Z = Zqd, Zqd/3, Zgd/10) suggest that such a diagram is not 
well suited to study and compare expansion properties of bright 
PNe in different environments since all tracks cover or embrace 
the observational data with no obvious trend with metal content. 
Moreover, at their maximum M(5007) brightness, all models, in- 
dependently of their metallicity (cf. Fig. 26 1, provide practically 
the same halfwidth velocity of Vhwhm - 18 kms~', in astonish- 
ingly good agreement with the observations! 

The theoretically expected distribution of the HWHM ve- 
locities depends on the evolutionary speed through the Vhwhm- 
M(5007) plane and the central-star mass distribution. In order 
to illustrate the former better, we show in Figs. [32] and [33] his- 
tograms which give the visibility time At spent by a particular 
model within a velocity bin of A Vhwhm = 1 km s _I . Thus, At is 
a proxy for the expected number of PNe per velocity interval of 
1 kms~'. These histograms combine the information contained 



in Figs. 26 and 31 



In particular, Fig. 32 displays the results for the 0.625 Mq 



sequences. The total simulation time shown in these figures is 
much longer than seen in Fig. 31 and the tracks extend thus 



to fainter magnitudes. The Vhwhm bins are additionally broken 
down into magnitude intervals to account for the fact that the 
velocity is not a monotonic function of M(5007) (Fig. 31 1. For 



instance, the low-velocity stage is passed rather quickly, which 



may explain the paucity of low-velocity objects seen in Fig. 31 



We show, for comparison, also the R out distributions. They 
appear bimodal where higher velocities belong to the high- 
luminosity phase as the central star crosses the HR diagram, and 
lower velocities to the recombination stage as the stellar lumi- 
nosity fades. 

The histograms seen in Fig. [33] belong to three simulations 
with the 0.696 Mq central-star model. Although the time scales 
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Fig. 32. Times At spent by the models within the velocity intervals A Vhwhm = 1 kms and AR out = 1 kms . The Vhwhm bins 
are broken down into M(5007) intervals as indicated by the different shades. The AR out histogram serves only for comparison and 
is not broken down into magnitude intervals. Shown are the three 0.625 M0 sequences with Z = Zgd (left), Z = Zgd/3 (middle), 
and Z = Zgd/ 10 (right) whose velocity and brightness evolution are also seen in the Figs. 26 and 31 The maximum post-AGB age 
considered is 5000 years. The At for the faintest magnitudes are lower limits only because of this time limit. 
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Fig. 33. The same as in Fig. 32 but for the three 0.696 M0 sequences with Z = Zgd (left), Z = Zgd/3 ( middle), and Z = Zgd/ 10 



(right) whose brightness evolution is also seen in Fig.|3T| The maximum post-AGB age considered is 1200 years. The Af for the 
faintest magnitudes are lower limits only because of this time limit. 



are shorter, the general behaviour is the same as already dis- 
cussed above. The histograms shown in Figs. 32 and ^demon- 



strate clearly that in a magnitude limited sample in which only 
the brightest PNe are contained all Vhwhm are quite low and 
cluster around 20 kms -1 . Very low halfwidth velocities (<10 
kms -1 ) occur only for faint models at the beginning of their 
brightness evolution (cf. Fig. [26}. Considering only models 
within the brightest bin of each sequence, the means of their 
halfwidth velocities are between 18 and 20 kms -1 , i.e. virtually 
independent of their metal content. 



From the discussion of the Figs. 31-33 we see that our mod 



els are in excellent agreement with measurements of bright PNe 
in distant stellar systems. The Vhwhm distribution of bright 
Galactic bulge PNe peaks between 15 and 20 km s -1 ( Gesicki & 
Zijlstra 2000 , Fig. 6 therein), and the 1 1 bright Virgo cluster PNe 
with spectroscopically resolved [O ra] 5007 A lines have Vhwhm 
between 12 and 22 kms -1 , with a mean value of 16.5 kms -1 
( |Arnaboldi eFaLl[2008l >. Th e bright PNe o f the Local Group of 
galaxies plotted in Fig. 4 of Richer ( 2006 ) cluster in the bright- 
ness range of M(5007) -4. 2 ... - 3.5, and their eye-estimated 
mean Vhwhm is - 18 kms -1 . Richer et al. ( 20 1 0b 1 argued again 



: 1 8 km s . Richer e 
that the mean 5007 A line width is not dependent upon maxi- 
mum [O ra] luminosity, metallicity or age of the progenitor popu- 
lation. All these findings are well explained by our hydrodynam- 
ical models, especially by the 0.625 Mq sequences (Fig. 31 1. 



All sequences discussed so far have the same initial model, 
viz. they start with an assumed constant AGB mass-loss rate of 
1 x 10 -4 M0yr -1 and a constant outflow velocity of 15 kms -1 . 
Here we briefly investigate how the above findings depend on 



initial mass-loss rate and wind velocity. For this purpose we 
computed additional sequences for the 0.696 M0 post-AGB stel- 
lar model (a = 2): The initial mass-loss rate of 1 x 10 -4 Mq yr -1 
was doubled/halved and the initial wind velocity of 15 krns -1 
lowered/increased by 5 km s -1 . The results are shown in Fig. 



34 



The interpretation of the influence of the AGB wind ve- 
locity is simple: the initial velocity enters as a corresponding 
up and downshift of Vhwhm(5007), rather independently of the 
metal content (Fig. 34 bottom, sequences 2/3 V ag b with Zgd and 



Zgd/3, and 4/3V ag b7~ZGD)- AGB mass-loss rate differences have 
a more subtle impact. The increased M ag b leads to a correspond- 
ing lower value of Vhwhm(5007) at the luminosity maximum, as 
expected. A lower M ag b, however, does not change the expansion 
much at maximum brightness (Fig. 34 top). The reason is that, 



because of the less dense AGB wind, the nebular shell becomes, 
in this particular case, optically thin well before its maximum 
possible [O in] line luminosity is reached. The [O in] brightness 
becomes temporarily lower until the shell recombines and be- 
comes again optically thick, and hence brighter. Concerning the 
expansion behaviour at maximum brightness, we see also that 
the metal content is virtually of no influence (Fig. [34j top, se- 
quences 2M agb , Z GD and 2M agb , Z GD /3). 

For completeness we add the discussion of two sequences 
with initial envelopes generated by radiation-hydrodynamics 
simulation along the upper AGB (sequences Nos. 6 and 6a of 
|Paper IV} . Both sequences have normal metallicities, i.e. Z = 
Zgd, and their Vhwhm evolution with M(5007) has already been 
displayed in Fig. 31 The corresponding histograms are seen in 
Fig.[35] 
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The behaviour of these sequences is somewhat different: At 
Vhwhm(5007) a 10. . . 11 kms _I the models become optically 
thin and reach then also their maximum [Om] brightness. The 
HWHM velocity further increases steadily during the optically 
thin part of evolution until recombination begins. The total range 
of Vhwhm(5007) within the brightest 3 magnitudes is for both 
sequences ^4 to ^25 kms~'. On average, the real expansion, as 
determined by R oul , is again considerably higher, viz. nearly a 
factor of two (Fig. 35 1. 

Since the [Om] brightness depends on the nebular metal 
(oxygen) content, it might be more useful to plot Vhwhm(5007) 
against the H/3 brightness in order to get rid of metallicity ef- 
fects in both the observation and theory. We expect that in a 
Vhwhm-M(H/3) plane the tracks of models which differ only in 
their metal content will be very similar, at least during the early 
evolution. Figure 36 confirms our conjecture: there is a rather 
tight correlation between M(H/3) and Vhwhm, and at maximum 
H/3 brightness (M(H/3) - -2) the HWHM velocities are ^ 16 ± 2 
kms 1 for 5007 A and ^20 + 2 for H/3. Larger differences be- 
tween the tracks occur during later phases when recombination 
determines the ionisation structure of our nebular models for 



ages above 3500 years (cf. Fig. 26 1 



The moderate variation of Vhwhm(HjS) is easily explained: 
The hydrogen line emerges from the entire ionised region, and 
its line width is controlled by the fast moving gas behind the 
outer shock (cf. Fig. [24] top panels) together with the (compara- 
tively) high thermal broadening (» 10 kms -1 ). In the beginning, 
thermal broadening dominates, but later expansion broadening 
becomes more important. Thermal broadening of oxygen lines 
is four times less, and the velocity evolution of, e.g., [O ra] 5007 
A can reasonably be traced by the line width. Very early, the 2+ 
zone is restricted to the innermost nebular regions only, close to 
the contact discontinuity where the gas is nearly stalling. Hence, 
Vhwhm(5007) is very low, much lower than the AGB wind ve- 
locity: only 4-5 kms~' as compared to V ag b = 15 kms _I . Later, 
Vhwhm(5007) steadily increases mainly because the 2+ zone 
extends until it embraces the whole nebula and approaches fi- 
nally V H WHM(Hy6). 

Our model simulations suggest that large differences be- 
tween velocities measured from the widths of hydrogen lines 
and lines of heavier ions are not only the result of different ther- 
mal broadening. Rather, these differences also point to the exis- 
tence of a substantial positive velocity gradient within all PNe. 
Such velocity gradients, as they are predicted by hydrodynami- 
cal simulations, have indeed been found for a few galactic PNe 
by means of a careful analysis of high-resolution slit spectro- 
grams (cf., e.g., |Gesicki etal.|1996|p58) . 

Recently, Rich er et al.| ([2009) presented halfwidth velocities 
of Galactic bulge PNe based on both Ha and [Om] 5007 A. 
These measurements are shown in Fig. 37 together with selected 



model sequences relevant for bright bulge objects of the Milky 
Way. One sees that the HWHM velocities predicted by our mod- 
els are in excellent agreement with measurements of the slowly 
expanding objects with unsplit lines, which are obviously also 
the younger ones of this sample (filled symbols). The agreement 
becomes poorer for a majority of the larger, spatially resolved 
PNe with split lines (open symbols). 



However, for a full appreciation of Fig. 37 we have to con 



sider that all Galactic bulge objects are to some extent spa- 
tially resolved. The typical coverage of the slit is between 30% 
and 50% for the slowly expanding objects (filled circles in 
Fig. 37 1 and usually less for the faster expanding PNe with dou- 
bly peaked profiles. Going from low to high spatial resolution 



increases the width of the line profiles of both the [O in] and H/3 
line by very similar amounts: In the low velocity regime we es- 
timated from model calculations typical values of 3 kms~' for 
5007 A and 4 km s _1 for H/3. For doubly peaked profiles this shift 
is larger, and this effect may explain the slight offset to the right 
that is seen at higher velocities between the observations and our 
models (open symbols in Fig. 37 1. 



Despite of the small inconsistency between the velocities de- 
rived from our (unresolved) theoretical line profiles and those 
measured by Richer et al. (2009 2010a), the interpretation of 



Fig. 37 is not really affected: We see how the ionisation stratifi- 
cation develops with time (or velocity) in nebular structures with 
positive velocity gradients in the same way as appropriate hy- 
drodynamical models predict. Note also that the two metal-poor 
sequences (Zgd/3, Zgd/10) do not predict the correct relation 
between hydrogen and oxygen line widths and are too faint in 
5007 A at low velocities. 

Spectroscopy done by |G6rny~ et al. ( 20 04|> reve aled that the 
sample of bulge PNe used by |Richer et al.| $M)\ |2010a| l con- 
tains several objects with WC central stars, which are marked 
in Fig. 37 by crossed symbols. Remarkably, their velocity be- 



haviour appears to be the same as that of their 'normal' counter- 
parts (apart from the outlier). Of course, the number of WC-type 
objects in this Richer sample is too small in order to allow a 
more definitive statement. 

Richer's linewidths measurements of PNe of the Galactic 



bulge are finally plotted in Fig. 38 over the H/3 absolute bright- 
nesses, assuming a distance to the bulge of 7.5 kpc and selecting 
only the younger objects with unsplit lines. There is again close 
agreement with the predictions of our model sequences. Central 
stars of about 0.6 Mq appear to be sufficient to explain the ob- 
served maximum H/3 brightness achieved by the PN population 
of the Galactic bulge. 

We can summarise the results of this section as follows: 

- Measuring the expansion of bright extragalactic PNe by 
means of the HWHM method, the only possible one for dis- 
tant objects which cannot be spatially resolved, is not very 
fruitful for more detailed studies. Although there exists an 
evolution of Vhwhm with time in, e.g. [Om] lines, the de- 
pendence on metallicity is extremely small and can virtu- 
ally be neglected. The time evolution of Vhwhm does by no 



means reflect the real expansion (cf. also Fig. 26 1. Typical 



HWHM velocities at maximum nebular brightness are SS20 
kms -1 , for all metallicities considered here. Lower speeds 
seen during the early PN evolution are entirely due to ioni- 
sation stratification in conjunction with a positive radial ve- 
locity gradient. However, the real expansion speed depends 
on metallicity and can be larger by factors of two or more 



(Figs. 26 32 and [33] see also the next section). 
Based on our hydrodynamical simulations reported here we 
state further that the very narrow range of halfwidth veloci- 
ties (Vhwhm) measured for the brightest PNe of galaxies with 
quite different properties as reported by Richer (2007 , Fig. 1 
therein) leads to only one conclusion: Stellar evolution at 
the tip of the AGB occurs independently of the properties 
of the underlying stellar population. The close agreement of 
observations with our model predictions implies that the ini- 
tial condition chosen by us at the tip of the AGB for, e.g., 
a 0.625 Mq stellar remnant, viz. M agb = 10~ 4 Moyr~' with 
V ag b = 15 kms -1 , are very close to reality. 
The narrow range of observed halfwidth velocities indicates 
further that everywhere the same main process for forma- 
tion, shaping, and evolution of PNe is at work: heating of 
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Fig. 39. flout /Vhwhm vs. Vhwhm 
for [Om] 5007 A {left) and 
H/3 (right) and for the differ- 
ent types of hydrodynamical se- 
quences. The evolution proceeds 
from low to high line widths, or 
Vhwhm- 



sufficiently dense circumstellar wind envelopes by photoion- 
isation, setting up an expanding shock wave which defines 
the outer edge of the PN. Wind interaction and/or magnetic 
fields, although important in single cases and/or during cer- 
tain phases, is obviously not the main shaping agent, con- 
trary to common wisdom. 



4.3. The real expansion 

Our simulations demonstrate clearly that it is impossible to de- 
rive the real expansion properties of PNe from a simple mea- 
surement of the halfwidths of spatially unresolved line profiles. 
Because of the ionisation and/or density structures together with 
projection effects the contribution of the fastest moving matter 
to (integrated) line profiles is strongly reduced. 

However, we can use our hydrodynamic models to quantify 
the relation between the real expansion, as defined by the prop- 
agation of the outer shock, and the spectroscopically measured 
expansion deduced from the HWHM of an emission line profile, 
which is in practice that of [O m] 5007 A, or also that of Ha 
or Hfj. Figure 39 illustrates how the ratio between shock speed 
and Vhwhm of [Om] 5007 A and H/3 (or Ho-) changes during 
the nebular evolution, as based on our simulations. We see that 
all types of models behave in the same way: The ratio between 
both velocities generally decreases with time and comes close to 
unity during the end of evolution. 

The behaviour of R out , 
plained as follows: 



can be ex- 



and Vhwhm seen in Fig. 39 i 
First of all, we recall that R out is determined 
by the density gradient of the (fully ionised) ambient medium 
and the sound speed of the nebular gas behind the shock front 



(cf. Paper III. Vhwhm, however, contains all the globally aver- 
aged information of gas motions and the density distribution in- 
side the nebula measurable via the Doppler effect. If the nebula 
has a stratified ionisation structure, Vhwhm depends also on the 
ion used. 

For instance, during the early ionising phase of the PN evo- 
lution the ratio between R oul and Vhwhm is quite large (greater 
than 3 for Vhwhm(5007) ^ 10 kms -1 ), since 2+ is restricted 
only to the strongly decelerated inner regions (cf. Fig. 
panel). The situation is different for hydrogen lines (Fig. 



13 



39 



right 
_ right 

panel): Since hydrogen ionisation starts while the matter is still 
moving with the AGB-wind velocity, Vhwhm starts with similar 
values of, e.g., 10-15 kms _I . The ratio of R out and Vhwhm does 
not exceed a factor of about two. 



With progressing ionisation, Vhwhm(5007) catches up later 
in the evolution because progressively faster moving matter from 
the outer regions behind the shock contributes to the integrated 
line profile. During the recombination phase it may happen that 
the shock is being decelerated so much that R out may fall even 
below Vhwhm, i.e. R m becomes lower than the flow speed of 
the nebular gas (cf. Fig. 39 left, the 0.696 Mq, Zqd models with 
V H whm(5007) > 30 km s^). 



We see especially in the left panel of Fig. 39 also that for 



fixed Vhwhm and for the same configurations (central star mass 
and initial envelope model) the correction for R out increases 
somewhat if the metallicity decreases. The shock propagation 
speed is controlled by the electron temperature and makes metal- 
poor and hotter nebula more dilute in their outer regions. The 
inner denser regions, however, are less accelerated because of 
the reduced wind power. Altogether, the width of the integrated 
line profile does not change to the same extent (full ionisation 
provided, cf. Fig. 26 1. - Such a metallicity influence on the line 
width is not seen for hydrogen lines (Fig. 39 right), most likely 
because of the generally larger thermal contribution to the line 
broadening. 

The interesting phase from an observational point of 
view occurs when nebulae evolve through their brightest 
stage (in [Om]). We have seen above that this happens for 
Vhwhm(5007) =! 15 ... 25 km s _1 , which implies correction fac- 
tors of 1.5.. . 2.5 to get the real expansion speeds. A mean cor- 
rection factor of two appears therefore to be a reasonable choice 



for the brightest nebulae of a sample. Figure 39 (right panel) 
suggests that two is also a quite useful mean correction factor 
for hydrogen lines. 

4.4. The case of the Magellanic Cloud PNe 

The Magellanic Clouds (MC) contain a rather large population 
of planetary nebulae which can also be studied down to fainter 
magnitudes than it is possible in most other extragalactic sys- 
tems. Thus, the MC sample of PNe is well suited for studying 
aspects of PN evolution in an extragalactic system. 

|Dopita et al.| ( |1985| |1988| l used in their seminal work 
of the Magellanic Cloud PNe a somewhat different defini- 
tion to measure the expansion of a spherical gaseous shell: 
V exp = 0.911 Vfwhm, where Vfwhm is twice the HWHM veloc- 
ity. Assuming a Gaussian shape for the intrinsic profiles, this is 
equivalent of V exp being half of the full width at 10% maximum 
line intensity, V hw %iom, or V exp = 1.82 Vhwhm, with Vhwhm be- 
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ing the velocity deduced from the HWHM of the line profile. 
Thus, Dopita et al.'s "expansion" velocities, V exp , are higher by 
nearly a factor of two as compared to the more widespread use 
of Vhwhm- The advantage of this definition is, however, that the 
such measured velocity is much closer to the maximum gas ve- 
locity behind the outer shock (see the discussion in Sect 

an d Table |3|>. 

Dopita et al. ( 1985 1988)) did not measure Vhw%iom directly, 



3.2.1 



but used Vexp = 1.82 Vhwhm instead, thereby implying implicitly 
a Gaussian shape for the intrinsic line profile. In reality, however, 
the spatially integrated line profiles are by no means Gauss func- 



tions (cf. Fig. 24 bottom panels). A moderate broadening by fi- 
nite spectral resolution does not change this. Additionally, the 
shape of the line profile changes with time because the density 
distribution and the velocity field develop as the model traverses 
the Hertzsprung-Russell diagram. 

The whole situation is illustrated in Fig. |40] where the rela- 
tions between R out and Vhwhm (top panel) and between Vhwio%m 
and 0.911Vfwhm (bottom panel) are displayed as a function of 
the nebular brightness for the different hydrodynamical PN se- 
quences used here. Vhwio%m expresses the velocity derived from 
the 10% level of the line profile. For a Gaussian we expect 
Vhwio%m/0.911Vfwhm = 1 (see above). From Fig. [40] (bottom 



panel) we see that the line profiles deviate considerably from a 
Gaussian shape, mainly during the late evolution when the mod- 
els become more extended. 



The top panel of Fig. 40 is equivalent to Fig. 39 but shows 
more clearly how the correction factor, to be applied to Vhwhm. 
depends on the nebular brightness. One sees again that a factor of 
two is quite reasonable for most models while evolving through 
their brightest stage. One sees also more clearly that only for 
the metal-rich model sequence around 0.696 Mq this factor falls 
below unity along the fading branch of evolution, as seen already 
above (Sect. |4~3] >. 

4.4.1. Nebular excitation and expansion 

From the rather large samples of LMC and SMC PNe|Dopita et| 



al. ( 1985 1 and Dopita et al. ( 1988 ) found a positive correlation 



of V exp with the nebular excitation parameter E, defined as 



E = 0.45 [X(5007)/X(H/?)] , 

E = 5.54 [X(4686)/X(H/3) + 0.78] , 



0.0<£<5.0, 
5.0 < E < 12.0, (5) 



where X stands for line luminosities or fluxes, respectively 
fDopita et al.|19 92). This correlation must be interpreted as the 
evolution of the nebular expansion speed with time or stellar ef- 
fective temperature. 

Equations Q constitute a purely empirical definition of 
the nebular excitation state in close resemblance to the 
original |Morgan| (|1984| l classification scheme (cf. Dopita & 
|Meath eringham]" fi990| ). Although the low excitation branch, 
E < 5, depends directly on the oxygen abundance, Eqs. <|3j are 
quite useful as a proxy for the nebular evolution. However, E 
is, like the brightness, neither a unique function of time nor of 
stellar temperature (see, e.g., Fig. 15 in Paper IV I. 

To perform a new interpretation by means of our hydrody- 
namic sequences, we proceeded as follows: Firstly, we computed 
the excitation parameters E of LMC and SMC PNe from the 
line fluxes available in the literature by means of Eq. (|5) and 



used V eX p as determined by Dopita et al. ( 1985 1 andlDopita et al 
1988). We discarded objects whose line profiles could only be 



models. We did then the same with our nebular models, i.e we 
derived E from the model's line strengths and determined V exp 
from the line profiles according to the Dop ita et al. [ prescription 
given above and applying an appropriate broadening (Gaussian 
with a FWHM of 10 kms~') to simulate instrumental broaden- 
ing. 



The result is shown in Fig. 41 which provides not only veloc- 
ity and brightness information of Dop ita et al.| ([T985 1988 ), but 
also oxygen abundances determined by photoionisation mod- 
elling of Dopita & Meatheringham ( 1991a|b[ l. The left panel of 
Fig. |4TJ displays the Magellanic Clouds objects, broken down 
into classes of different M(5007) brightnesses and oxygen abun- 
dances (relative to hydrogen), as indicated by the symbol sizes 
and their grey values. This panel is a new version of Fig. 4 in 
Dopita et al. (1988 ) with more accurate excitation values. 



There is a broad positive correlation between excitation class 
and expansion velocity of the Magellanic Clouds objects in the 
sense that there are no low velocities for high excitation and 
no high velocities for very low excitation. The distribution is 
very narrow at low velocities (V exp S5 15 kms~'), but becomes 
rather wide at higher velocities or excitations, independently of 
the oxygen abundance. One notices, however, that the brightest 
objects have preferentially also the highest amount of oxygen. 



The left panel of Fig. 41 shows again what is already evident 
from our Fig.[3T]and Fig. 4 in |Dopita et al.| ( [T992| l: The most lu- 
minous objects, which accumulate around E ^ 6 (big symbols in 
Fig. 41 1, have rather modest expansion rates (in terms of Vhwhm, 
top abscissa) between 15 and 20 kms~'. The fainter objects 
have, in general, either low (young objects) or high expansion 
rates (older objects). Obviously PNe in the Magellanic Clouds 
show the same expansion behaviour as PNe in the Galactic 
bulge, the Local Group, and the Virgo cluster. 



The middle panel of Fig. 41 compares the observations with 



the predictions of our hydrodynamical simulations in terms of 
excitation, brightness, and halfwidth velocities. We used the 
slightly metal-deficient models with Z = Zgd/3 in order to take 
the metallicities of both Clouds roughly into account. Note that 
the excitation E, if computed according to Eqs. Q, does not 
reach the limit E — 5 from below because of the reduced oxygen 
content, creating thereby an artificial gap in the presentations of 

Fig- ED 

For a full appreciation of this figure it is important to realise 
that the excitation E does not only depend on the stellar tem- 
perature but also on the nebular density. Lower densities favour 
higher excitations for a given stellar temperature. Therefore, the 
models around the 0.605 Mq central star, which are less dense 
than the models around the 0.696 Mq star and become optically 
thin during the course of evolution, reach higher excitations than 
the 0.696 Mq models In spite of the fact that for the latter the 
central star achieves much higher temperatures (e.g. 240000 K 
instead of only 160 000 K for the 0.605 Mq case). We showed in 
Paper IV (Figs. 13 and 14 therein) that the highest degree of ex- 
citation (E ss 10) is only reached by optically very thin nebulae. 



fitted by a double Gaussian since they have most likely an ex- 
treme bipolar shape not suitable for comparison with spherical 



An example is the a = 3, 0.595 Mq sequence shown in Fig. 41 
(middle): It does not recombine much and reaches the highest 
expansion rates at the highest possible excitation level, as given 
by the He/H ratio which is 0. 1 1 for our models. 

Figure [4T| demonstrates that all objects start their expansion 
(at E a 0), as the models do, very slowly, with Vhwhm(5007) of 
only a few km s _1 . Since our models have AGB-wind flows with 
10 or 15 km s _! , these low values are thus no indication for a sim- 
ilarity low AGB-wind velocity. We have already explained above 
this phenomenon as the consequence of the ionisation stratifica- 
tion (of oxygen) in combination with a steep velocity gradient 
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Fig. 41. Nebular excitations, E, vs. 5007 A halfwidth velocities, V exp = 0.9HVfwhm, for Magellanic Clouds PNe (Dopita et al. 
19851 1988] > and our hydrodynamical simulations. Observed objects and models are broken down into the same M(5007) bright- 
ness classes, as indicated by the symbol sizes: the larger the symbol, the brighter the object or model, respectively. The approxi- 
mate brightness range covers 4 magnitudes. The top abscissa is labelled in terms of Vhwhm for easier comparison with previous 
figures. Left panel: LMC and SMC PNe additionally broken down into classes of different oxygen abundances (from Dopita & 
Meatheringham| 199 1 a|b| l, visualised by various grey shades. The darkest grey corresponds to >Zgd/2, the lightest grey to <Zgd/5. 
The vertical arrow indicates the oxygen abundance corresponding to Zgd/3. Middle panel: comparison of observations with the 
model sequences as indicated. The gaps seen in the model sequences are artefacts due to the definition of E combined with the 
lower oxygen abundance of the Clouds (Eqs.|5] see text there). Note that the spacings of the models within a sequence does by no 
means reflect the speed of evolution. Right panel: A 0.625 M0 model sequence (a = 2, Z = Zgd/3) with different initial parameters, 
M dg b = 0.5 x 10~ 4 Moyr -1 and V ag h = 7.5 kmr 1 (filled triangles), compared with the 0.625 Mq Zgd/3 sequence from the middle 
panel and a metal richer (Z = Zgd) 0.625 M0 sequence (see legend). 
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within the ionised shell during the early evolution of a PN (cf. 
Fig.fTJl. 



During the course of evolution the line widths increase, with 
some spread in terms of excitation and Vhwhm which is caused 
by the different combinations of central star and circumstellar 
envelope. Obviously, Fig. HT] (middle) implies again that for all 



observed PNe the same mechanism for their formation and ex- 
pansion must be at work, which, according to our models, is 
heating by ionisation and the immediate creation of an expand- 
ing shell bounded by a shock front. 

During recombination, E drops to intermediate values al- 
though the star remains very hot (T e s ^ 120000 K). The ioni- 
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sation becomes stratified, and the 2+ is confined to the inner re- 
gions closer to the hot central star. The outer shock is decelerated 
because of the post-shock pressure drop due to recombination 
while the inner region is still accelerated by the thermal pressure 
of the hot bubble (shocked stellar wind gas) and somewhat later, 
by reionisation. Because of this, Vhwhm may decreases during 
recombination, but continues later to increase further (cf. 0.605 
Mq sequence in Fig. 41 middle panel). 

From Fig. 41 (middle panel) we see, however, that our model 
sequences do not match the observations perfectly, which, of 
course, is of no surprise since the models are not tuned to match 
objects of the Magellanic Clouds. For the most relevant se- 
quence, 0.625 M0 and Z = Zgd/3, E falls behind the Vhwhm 
evolution, while the 0.605 Mq sequence has the correct expan- 
sion property but remains a little bit too faint at medium exci- 
tations. The 0.595 Mq sequence stays always far too faint (cf. 
their small symbols) and is thus omitted for further considera- 
tions. The 0.696 Mq sequence becomes somewhat too bright, 
although the Vhwhm evolution is marginally right. 

We have already seen in Fig. [34] that a variation of the ini- 
tially chosen AGB-wind velocity leads to a corresponding varia- 
tion of the PN expansion rate. Hence, we expect that models with 
different AGB-wind velocities would be shifted horizontally in a 
diagram such as displayed in Fig. 41 To test this possibility, we 
used the 0.625 Mq, Z = Z CD /3 sequence with M agb = 0.5 x 10~ 4 
Moyr -1 and V ag b = 7.5 kms~'. This combination of mass-loss 
rate and AGB-wind velocity ensures that the initial circumstel- 
lar density profile is the same as in the comparison sequence. 
However, because of the more slowly moving AGB material, 
also Vhwhm increases now more in pace with E, although ob 
jects with the highest excitations are now not covered (Fig. 41 
right panel). 

For comparison, we replotted here also the 0.625 Mq se- 
quence from the middle panel, supplemented by a corresponding 
metal richer sequence with Z = Zgd- One sees that just changing 
the metallicity will not lead necessarily to a better match with the 
observations. 

One concludes from this exercise that a mixture of models 
with somewhat lower final mass-loss rates and lower AGB wind 
speeds and of models as they are used here would lead to a good 
match with the observed HWHM velocity-excitation correlation 
of PNe in the MCs. 



4.4.2. Stellar temperature and nebular expansion 

Because of the dependence of E on both stellar temperature 
and (mean) nebular density, it is not surprising that the ob- 
served sample cannot be described by a single evolutionary 
track. A more direct comparison can be made if we correlate 
V ex p (= 0.991 Vhwhm) of the nebula with the effective tempera- 



ture of the central star (cf. Fig. 6 in Dopita & Meafheringham 
1 99 1 a) > . We used the same sample shown in Fig.|4T|and selected 
all those objects for which Dopita & Meatheringham | ( [l 99 1 a|b| l 
provided stellar temperatures, based on detailed photoionisation 
models. 

The use of the stellar temperature (if accurately known) as a 
proxy for the time evolution across the Hertzsprung-Russell dia- 
gram is useful since (i) it is a distance-independent quantity, and 
(ii) its range depends weakly on the stellar mass. Furthermore, 
the two driving forces for the nebulae, the radiation field and the 
wind, depend directly on the stellar effective temperature. 



that the signatures of nebular expansion, as measured by a suit- 
able line width, evolve with time (or stellar effective tempera- 
ture), and that this time evolution is correctly described by our 
hydrodynamical model sequences. Note that the theoretical se- 
quences shown in Fig.[42]cover very different time spans, rang- 
ing from about 10000 years for 0.605 Mq, through about 4000 
years for 0.625 Mq, to about 1000 years only for 0.696 Mq. 
Yet the sequences are nearly 'degenerate' if the stellar tempera- 
ture is used as a proxy for the evolutionary progress across the 
Hertzsprung-Russell diagram. 

A similar correlation between nebular expansion and stellar 
temperature was also found for Galactic PNe with double shell 
structures in our Pa per II| where typical flow speeds of shell and 
rim gas were measured from spatially resolved high-resolution 
echelle spectra. 

4.4.3. Nebular ionisation stratification and expansion 

An indication of the ionisation structure and the velocity field of 
a PN can be gained by measuring the HWHM velocities of dif- 
ferent ions of the same element, viz. of + and 2+ (e.g. Dopita 
et al. 1988, their Fig. 3). These authors found that velocities de- 
rived using + are in most cases higher than those from 2+ . We 
replotted their data in Fig. [42] (middle), supplemented again by 
our dynamical models in order to check whether our models are 
able to predict the same behaviour. 



We see from Fig. 42 (middle) that the observed trend of 
the expansion signatures are, indeed, quite well matched by 
our models, although the [Ora] velocities of the younger 
models are somewhat lower than observed. At the begin- 
ning of the evolution all models behave similarily, with 
Vhwhm([Oii]) > Vhwhm([0 in]), because of the radial ionisa- 
tion stratification of oxygen and the positive velocity gradient. 
Later, if the models are optically thin and 2+ becomes the main 
ionisation stage throughout the model, both ions, + and 2+ , 
trace about the same nebular regions, albeit with different frac- 
tions, and hence V H whm([Oii]) - V H whm([0 in]). This transi- 
tion occurs quite rapidly and produces the left turn of the model 



tracks seen in Fig. 42 (middle). 

Note that the models around the 0.696 Mq central star 
remain optically thick throughout the whole evolution, with 
Vhwhm([0 ii]) always greater than Vhwhm([0 hi]) until the star 
begins to fade. This sequence with its quite massive central 
star model, 0.696 Mq, is obviously not an option for PNe in 
the Magellanic Clouds, also because, as we mentioned already 
above, the models become too bright. 



Figure 42 (left panel) is an adaption of the Dopita & 
Meatheringham figure mentioned above and demonstrates again 



The interpretation of Fig. 42 (middle) is the following: If 
Vhwhm([Oii] is significantly greater than Vhwhm([0 in]), the 
object in question is ionisation stratified and optically thick, ir- 
respective of the amount of Vhwhm - If both velocities are similar, 
the object is optically thin, which is the case for most of the MC 
objects plotted in Fig. |42] (middle). It appears, however, that there 
exists often a tendency for Vhwhm([0 ii]) of being greater than 
Vhwhm([0 in]) by « 10 kms~', on the average. This behaviour 
is best matched by our 0.605 Mq sequence during the optically 
thin phase. 

4.4.4. Evolution of nebular masses 

As PNe expand into the AGB matter, their mass must increase 
with time. This has been convincingly demonstrated by Dopita 
|& Me atheringham ( 1991a|b i in their study of the Magellanic 
Cloud PNe. The mass increase with time (or radius) is ruled by 
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the speed of the (outer) shock and the density profile of the for- 
mer AGB wind matter. A decrease of the ionised mass is possible 
during the recombination stage as the central star fades. 

In the right panel of Fig. |42] we compare the |Dopita &| 
Meatheringham ( |1991a|b l data with our model predictions. We 



its maximum H/3 brightness well before the m aximum in 
5007 A is reached (compare Figs. 9 and 1 1 in Paper IV i, 



used the reduced mass as defined by Dopita & Meatheringham 



Before interpreting this figure one should say that nebular 
masses are difficult to determine, even if the distances are well 
known. Uncertainties of electron densities and radii can add up 
to quite large mass errors. Having this in mind one can only state 
that our simulations provide a rather satisfactory agreement with 
the observations. 



4.4.5. Nebular brightness and expansion 

We conclude our treatment of Magellanic Cloud PNe by ap- 
plying VHWHM-magnitude diagrams already introduced and dis- 

We thus coll ected the velocity infor- 
' Dopita et aT] ( |1985| l and |Dopita et al. 



cussed before in Sect. 4.2 
mation on MC PNe from 

(1988), converted their 'expansion velocity' values, V exp , in or 



der to be consistent with, e.g., Figs. 31 and 36 to the ones used 
here (Vhwhm = V exp / 1.822), and plotted V H whm(5007) against 
M(5007) and M(H(3), respectively (Fig. 03}. The model se- 
quences selected are the same as used in the previous figures, 
supplemented by sequences with lower (Zgd/10) and higher 
metallicity (Zgd)- Note that no unique relation exists between 
M(5007) and M(H/3): Each sequence has its own relation, de- 
pending on central star mass and corresponding nebular ionisa- 
tion. 

The left panel of Fig. 43 is the same as Fig. Fig. 



31 



but 

only for Magellanic Clouds objects, the right panel is a virtu- 
ally metal-independent presentation. This is seen by the close 
similarity between the Z = Zgd, Zgd/3, and Zgd/10 0.625 Mq 
sequences. Again, the observed distributions of the MC objects 
are very well matched by our model sequences. Several facts are 
evident from Fig. 43 



- The brightest objects are obviously explained by models 
with more massive central stars (^0.6-0.7 M0) and optically 
thick (or partly thick) nebulae, while objects with medium 
brightnesses need optically thin models with less massive 
central stars ($0.6 M0). 

- The M(5007) peak brightness of the SMC sample appears to 
be fainter than that of the LMC by about 0.2 mag, obviously 
because of the lower LMC metallicity (more precisely, oxy- 
gen abundance). This difference is not seen for M(H/3), as 
one would expect. The bright cut-off occurs at M(Hf3) -2 
mag for both clouds. 

- The model sequence with V ag b = 7.5 kms -1 can explain the 
bright objects with low V H whm(5007), but fails to explain 
those PNe with V H whm(5007) ^ 20 kms 1 . 

- The_jnean_VHWHM(5007) velocities for all objects of the 
Dopita et al. samples are ^17 km s (LMC) and ^15 km s _I 



(SMC). Considering that the SMC sample is rather small and 
that the (formal) standard deviations are about 6 kms~', we 
conclude that both mean values are virtually equal within the 
errors, in good agreement with our simulations. 
- Looking only at objects within the brightest magnitude in- 
terval, those which are the brightest ones in H/3 have, on the 
average, a lower halfwidth velocity than the corresponding 
[Om] bright objects: The mean Vhwhm (5007) for the for- 
mer is about 12 kms~' (right panel of Fig. 43 1, while the 
mean for the latter is about 18 kms~' (left panel). The rea- 
son for this difference is simply the fact that (i) a PN reaches 



at a phase where the expansion is still rather slow, and that 
(ii) most PNe become optically thin and thus fainter in H/3 
shortly afterwards. We remind the reader that the nebular ex- 
pansion rate generally increases with time, especially during 
the early, optically thick stage. 



5. Discussion 

The work presented here is the first where modern radiation- 
hydrodynamics simulations of the evolution of planetary neb- 
ulae are systematically applied for objects in distant stellar 
populations. To this end, we present a set of lD-radiation- 
hydrodynamics simulations which are especially tailored to 
study the formation and evolution of planetary nebulae in galax- 
ies and to interprete the observations. The model sequences 
cover a wide range of metallicity in terms of abundance distribu- 
tions scaled with respect to that which is typical for the Galactic 
disk. The different abundance patterns generally found within a 
galaxy are thus neglected here. Important metal-dependent pro- 
cesses (e.g. line cooling) are considered in detail, but the depen- 
dence of the central-star wind on metallicity, which is of rele- 
vance for the PN shaping process, could only be treated approx- 
imately. Direct consequences of the stellar chemistry on details 
of the post- AGB evolution are not considered. 

The initial conditions in terms of AGB mass-loss rate and 
wind speed, and the range of masses considered for the cen- 
tral stars, are chosen such that (i) the emerging model structures 
mimic closely the observations known for Galactic disk objects, 
and that (ii) the bright end of the PN luminosity function can 
be explained by optically thick model nebulae around central- 
star models with ^0.63 M0. Typical values needed for the final 
AGB mass-loss rates are M ag b - 10~ 4 Mq yr~' with wind veloc- 
ities between about 10-15 km s _I . Values of about 10~ 4 Mq yr~' 
for the final AGB mass-loss rates from four well-known Galactic 



PNe have recently been derived by Sandin et al. ( 2008| l using in- 
tegral field spectroscopy of the respective PN haloes. 

In order to ease the comparison of our models with dis- 
tant PNe which are usually (except for the Galactic bulge) not 
spatially resolved, observable quantities are integrated over the 
whole model nebula. Emphasis is given to the extraction of ve- 
locity information from the half widths of emission line pro- 
files (Vhwhm)- Morisset & Stasiriska (2008) showed recently by 
means of non-spherical model nebulae that HWHM-velocities 
derived from the profiles of volume integrated emission lines are 
expected to be a robust measurement in the sense that these ve- 
locities do not depend much on morphology and orientation of 
the model, and also not much on the details of the velocity law. 
Spherical models as they are used here appear thus fully justified 
for deriving mean properties of a larger sample of objects. 

We emphasise that all existing efforts to model axisymmet- 
ric PN structures by fast collimated winds (jets) and/or mag- 
netic fields are using initial AGB wind configurations that are 
far from reality, at least for the typical PNe discussed here. The 
initial AGB-wind densities employed are always very low (typ- 
ical values correspond to M ag b ~ 10~ 6 Mgyr -1 ) so that ther- 
mal pressure is dynamically unimportant (see. 



e.g.. 



Akashi & 



|Soker|2008] l. However, as we have mentioned above, typical PN 
structures require final AGB-wind mass-loss rates of about 10 
Mq yr _1 , higher by about two orders of magnitude than assumed 
in current jet modelling. Additionally, the inner boundary condi- 
tions are fixed, i.e. no allowance for an evolving star and wind is 
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Fig. 43. Comparison of HWHM 
velocities vs. absolute magni- 
tudes of Magellanic Cloud PNe 
with different model predictions 
as indicated in the legends. The 
dashed-triple-dotted line corre- 
spond to the V ag b = 7.5 kms~' 
sequence. The velocity measure- 
ments are again from |Dopita et] 
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is shown, instead of their V exp , 
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made, although the simulation times often become comparable 
to the star's transition time from the tip of the AGB towards the 
PN region (« 10 3 years). 

We close the discussion with the warning that our simu- 
lations are exclusively based on the properties of hydrogen- 
burning central-star models and should not be used to explain the 
evolution of PNe around hydrogen-deficient (i.e. Wolf-Rayet) 
central stars. The samples used here for comparisons contain cer- 
tainly an unknown fraction of such objects because it is not pos- 
sible or feasible to sort them out. The mere fact that they do not 
seem to disturb the correlations found in this work may indicate 
that the fraction of PNe with Wolf-Rayet central stars is rather 
small and/or that their evolution is not much different from that 
of their counterparts with a hydrogen-rich central star. There ex- 
ists also the possibility that a different evolution of Wolf-Rayet 
central stars is hidden because of inappropriately chosen observ- 
ables such as Vhwhm- 

6. Conclusions 

Our extensive radiation-hydrodynamics simulations presented 
in this work show clearly that, contrary to common belief, 
wind interaction is not the main driving force for the formation 
and expansion of the bulk matter in planetary nebulae around 
hydrogen-burning central stars: The metal-poorest models with 
the weakest stellar wind expand fastest. Instead it is the pressure 
gradient built up by ionisation and heating of the circumstellar 
gas that drives a D-type shock through the ambient medium and 
creates thereby an expanding circumstellar structure consisting 
of different shells, called a PN. So to speak, it is the radiation 
field of the star, not its wind, that is the main agent for formation 
and evolution of a PN. The overall shape of a PN is then con- 
trolled by the matter distribution set up by the wind before ioni- 
sation starts. Possible processes have been proposed by Dijkstra 



& Speck (2006). The central-star wind is only necessary to pre- 
vent the circumstellar matter from falling back to the stellar sur- 
face and may compress the inner nebula to a denser shell, called 
the rim, when the wind power becomes strong enough at higher 
metallicities. 

Although these wind compressed rims are often the most 
conspicuous parts of a PN, they contain only a very small frac- 
tion of the total nebular mass. The outer boundary of the rim, 
i.e. the rim's shock, separates then the bulk of nebular matter 
which is thermally expanding. Using a stellar wind under con- 
ditions typical for Galactic disk objects, the rim gains no more 



than about 20% of the total ionised nebular matter at the end of 
the horizontal evolution. Only once recombination in the shell 
has started, the (ionised) mass contribution of the rim becomes 
much higher and may even become dominant. At low metallic- 
ities, the rim becomes unimportant at all and loses its typical 
signature. 

Based on our extensive modelling and existing observations 
the basic properties of PNe in different environments can be 
summarised as follows: 

- The metal content has a profound influence on the expansion 
properties and on the development of the overall morphol- 
ogy. For low metallicities, the expansion is fast due to the 
high electron temperatures while the wind interaction gen- 
erated from the central star is weak, leading to an only rel- 
atively low bubble pressure; hence low-metallicity PNe are 
expected to be rather extended. On the other hand, at high 
metallicities, expansion is slow and wind pressure strong; 
hence we see a much (geometrically) thinner structure with 
a prominent rim which is clearly separated from the outer 
shell matter. According to our models, the transition where 
the stellar wind becomes less important for shaping the in- 
ner parts of a PN occurs at about or slightly below Zgd/10, 
depending on the initial parameters of the models. 

- We showed that low-metallicity PNe are prone to devia- 
tions from thermal equilibrium because of their reduced 
line-cooling efficiency and fast expansion. This may already 
occur at about 1/10 of the solar (or Galactic disk) abun- 
dance distribution. Studies employing classical photoionisa- 
tion codes may become problematic in these cases. 

- In contrast to thermal equilibrium, the ionisation/recombi- 
nation equilibrium is always rather well established, except 
possibly for a very brief phase during the rapid stellar fad- 
ing after the 'horizontal' evolution through the HR diagram. 
Independent of metallicity, the halo of a PN is generally far 
from thermal and ionisation equilibrium. 

- Although the influence of metallicity on expansion property 
and density distribution is considerable, the half widths of 
(volume) integrated emission line profiles are rather insensi- 
tive to the metal content. The increase of the HWHM veloc- 
ity, Vhwhm (in, e.g., the 5007 A [Om] line) is initially due to 
the growing 2+ zone encompassing faster moving gas from 
the outer nebular regions, but later due to the generally in- 
creasing expansion speed as a consequence of sound speed 
behaviour and upstream radial density gradient. The evolu- 
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tion with time, or brightness, is such that all optically thick 
(or nearly optically thick) models used in this study reach, 
irrespective of their metallicity, a HWHM velocity of about 
18 kms -1 at their maximum [Om] brightness. Up and down 
shifts of this value are possible by changing the original AGB 
wind speed accordingly. 

Our model predictions are in remarkable agreement with re- 
cent observations presented by Richer et al. ( 201 Ob} which 
show a very limited range of the mean HWHM velocity, 
^15-20 kms , for a sample of the brightest PNe drawn 
from galaxies of the Local Group with quite different ages 
and metallicities. According to our modelling, these bright- 
est, optically thick PNe should have typical central-star 
masses of ^0.65 Mq, with final AGB mass-loss rates of 
^10 -4 M o yr -1 . 

The few investigated bright PNe in the Virgo cluster show 
the same behaviour as those in the Local Group: their mean 
HWHM expansion (in 5007 A) is ^16.5 kms -1 , implying 
that also the same parameters for the central stars and the 
final mass-loss rates hold for the bright PNe in the Virgo 
cluster. 

- It could be shown by our simulations that a more sensitive 
method to determine the PN kinematics must rely on mea- 
suring the expansion from the outermost wings of the line 

[T988> . 



profile, as already proposed by |Dopita et al. (1985 



Only in this case it is possible to get information about the 
fastest moving parts of a PN, usually those immediately be- 
hind the shock whose speed is also much more sensitive to 
the metal content than the HWHM velocity. This method re- 
quires, of course, high-quality line profiles which are still 
difficult to get for more distant stellar systems. 
- The run of nebular expansion with evolution (i.e. with time) 
cannot be determined by correlating a velocity measurement 
just with nebular (absolute) magnitude. One needs instead a 
parameter which depends directly on the evolutionary status 
of the central star, e.g. r en - . If T e g is not available, or too dif- 
ficult to determine, the nebular excitation, E, can be used as 
an distant independent proxy of the stellar effective temper- 
ature. 

The Vhwhm evolution seen for the PNe in the Magellanic 
Clouds, if plotted over E or T e g, are fully explained by our 
model sequences with the appropriate metallicity and with 
central star masses between =s0.6 . . . 0.7 Mq: Vhwhm(5007) 
ranges in both clouds from about 5 to about 30 kms -1 , but 
the brightest PNe (in 5007 AVhave HWHM velocities be- 
tween 15 and 20 kms -1 onlyp] Close to maximum H/3 the 
range is wider, ^5 . . . -20 km s , with a mean value of about 
12 kms -1 , again in full agreement with our model predic- 
tions. 

A similar velocity evolution with stellar effective tempera- 
ture (i.e. with time) was also found for Galactic disk objects; 
see Paper II| Fig. 12 therein, and also |Steffen & Schonberner| 



(2006). A direct comparison, however, is impossible because 
the velocities measured are based on spatially resolved line 
profiles and refer to the outer parts of the shell. They are usu- 
ally greater than the HWHM velocities of volume-integrated 
line profiles. 



Also |Richer et al. (T2008 ) found recently the signature of a 



modest velocity evolution by studying the Milky Way bulge 
PNe: the median line width in [Ora] 5007 A is 17.2 kms -1 
for the cooler central stars, but 23.7 kms -1 for the hotter, 



12 The Magellanic Clouds are included in the 
sample of galaxies. 



Richer et al. 



(2010b I 



more evolved ones. The standard deviations of both distribu- 
tions are 5.1 and 7.6 kms -1 , respectively. Note however that 
most of the bulge PNe are to some extent spatially resolved. 

- Another important result is that the HWHM velocities al- 
ways underestimate the true expansion speed of a planetary 
nebula. The amount depends on the evolutionary stage, but 
for objects near their maximum [Om] brightness, the real 
speed is about a factor of two higher (cf. Fig. 39 1. 

- The nearly identical expansion behaviour, (in terms of the 
line width Vhwhm) observed for the bright PNe in stellar 
populations which differ in age and metallicity allows the 
conclusion that most likely only one basic (and rather sim- 
ple) physical mechanism is responsible for driving the for- 
mation and expansion! According to our models this is al- 
ways photoionisation by the stellar radiation field. The stel- 
lar wind, as it is modelled here, does not play a dominant 
role. Also, it does not appear very likely that more 'exotic' 
scenarios like disk accretion and jets, or binary and common 
envelope evolution which are heavily debated, are of major 
importance, at least for the bright PNe with more massive 
central stars and envelopes. 

It has already be mentioned in the discussion above that 
our models of planetary nebulae are entirely based on radiation- 
hydrodynamics simulations of wind envelopes around hydrogen- 
burning central-star models. Some of these conclusions may not 
apply for planetary nebulae around hydrogen-poor, i.e. Wolf- 
Rayet, central stars. For instance, the winds of Wolf-Rayet cen- 
tral stars are much more powerful which makes them very im- 
portant for the formation and evolution of PNe. However, since 
the formation and evolution of hydrogen-poor central stars is still 
unclear, any detailed simulations of the sort as reported here are 
premature. 
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Fig. 20. Radial ionisation structures for hydrogen, helium, car- 
bon, nitrogen, oxygen, and neon of the Z = Zgd/10 model pre- 
sented in Fig. [7] The two top panels show again the underlying 
nebular structures: densities (ions = thick, electrons = dotted), 
gas velocity (thin), both left, and electron temperatures (dynam- 
ical = solid, equilibrium = dotted, right). The shadowed vertical 
strip marks the region occupied by the shell's leading shock in 
all panels. The ionisation fractions, X, are labelled as follows: 
dotted = neutral (hydrogen and helium only), short-dashed = 1st 
ionisation, dash-dotted = 2nd ionisation, dash-dot-dot-dot = 3nd 
ionisation, long-dashed = 4th ionisation, and dotted = 5th ioni- 
sation (except for hydrogen and helium), and are displayed in the 
left panels for the dynamical model (Xd yn ). The ionisation frac- 
tions are omitted inside the contact discontinuity (r < 3.5 x 10 17 
cm) in order to avoid confusion. On the right, the ratios between 
the equilibrium (eq) and dynamical case (dyn) are shown. Only 
the ratios of major ionisation fractions, which are seen also in 
the left panels, are rendered for clarity. 
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Fig. 21. The effect of dynamical/expansion cooling as described 
by the source/sink term p(V ■ v) for the Zgd/10, a — 3 model 
shown in Fig. [7] where the central star has gained its maximum 
effective temperature (see the track in Fig. H). Top row. model 
structure with particle densities and flow velocities (left, ions = 
thick solid, electrons = dotted, velocity = thin solid) and tem- 
peratures (right, dynamical = solid, equilibrium = dotted). The 
shadowed vertical strip indicates the thickness of the outer shock 
region. Second row: different contributions of line cooling (left) 
and dynamical cooling (right), as indicated in the legends, both 
based on the equilibrium temperature. The right panel with the 
two hydrodynamical terms (dv/dr and 2v/r) is "mirrored" at the 
horizontal axis in order to show also heating due to compres- 
sions if dv/dr becomes negative (shock front and halo/upstream 
region). Note that no cooling terms are plotted for the hot 
bubble, i.e. inside the contact discontinuity at ^3.5 x 10 17 cm. 
Bottom left (shadowed): total cooling (line + hydro) to line cool- 
ing, taken from the middle panels (thick solid) and for models 
with the same evolutionary time but with other metallicities (see 
legend). Note the different radius ranges covered by the models: 
the inner boundaries are set by the respective contact discon- 
tinuities, and the vertical 'double' lines mark the outer shock 
regions where heating by compression prevails. Bottom right: 
ratio between adiabatic cooling using Eq. Q, with the model 
data displayed in the top panels (hydro) and using various ap- 
proximations (approx), always based on the equilibrium elec- 
tron temperature: mean of total particle density of 40 cirT 3 and 
fixed outer radius R out = 1.4 x 10 18 cm (dash-dotted), the (same) 
mean density but a variable radius (dotted), variable density and 
fixed outer radius (dashed), and finally variable density and ra- 
dial coordinate (solid). In all cases, the model velocity Vhwhm = 
40 kms~' is used. 
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Fig. 23. Mean nebular [Ora] electron temperatures vs. nebular 
line ratio L(4686)/L(H/J) for the dynamical a = 3 sequences 
from Fig. 18 (dotted in thermal equilibrium). Again, the se- 



quences are only plotted until maximum stellar temperature in 
order to avoid confusion. The circles give the corresponding po- 
sitions of the Halo PNe: filled circles are based on the observed 
[Om] lines, open circles are results of detailed analyses using 
photoionisation models ( Howard et al.|1997] >. The inset renders 
an enlarged view of the region around L(4686)/L(Hy6) ^ 1. 
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Fig. 25. Evolution of nebular mass fractions vs. time for the 
0.625 M0 sequence from Fig. 24 with Z = Zqd- Shown are 



the total (ionised) mass enclosed by the outer nebular bound- 
ary, Mi on (R out ), with the different contributions of the shell (light 
grey) and rim (dark grey) indicated. The different lines indi- 
cate mass fractions contained within certain distances R from 
the star for which the real gas velocities reach distinct limits: 
M ion (R 



up 

HWHM 



) contains all mass shells with v < V^ HM (thick 



dashed line), and Mi on (R^^ M ) is the corresponding mass value 
for [Om] 5007 A (thick solid line). HWHM refers to the total, 
volume integrated line profile. The thin lines (solid and dashed) 
indicate the masses enclosed by the radii /?hwio%m where the 
(real) gas velocity v corresponds to the half width of the 10% 
intensity level of the integrated line profile, i.e. v = Vhwio%m- 
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Fig. 26. Three different expansion velocities of a PN: Vhwhm of 
the spatially integrated [O m] 5007 A line profile, the post-shock 
velocity, V p0 st, behind the outer shock at R out , and the shock ve- 
locity, R out , itself. Shown is the complete time evolution of these 
velocities for three sequences with the 0.625 M0 central-star 
model and different metallicities (see figure legend). The filled 
circles mark the snapshots shown in Fig. [24] which correspond 
also to the moment of maximum [O in] emission. The relevant 
model parameters are listed in Table [3] where also the values of 
Vhwhm, V post , and R out for the 3Zgd sequence, not plotted for 
clarity, are given. The bumps at age 1000 years are due to the 
optically thick/thin transition. 
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Fig. 27. The maximum efficiency of converting stellar radiation 
into [Oih] 5007 A line radiation, (L 50 07/i s tar) max , of 0.696 M Q 
(filled circles) and 0.625 M0 models (diamonds) against their 
metallicity, Z/Zq D . The upper abscissa gives the corresponding 
oxygen content relative to hydrogen, e(O) (= log(0/H) + 12). 
The numbers along the dotted line segments connecting the 
0.696 Mq data points indicate the slope. The inverted parabola 
(dashed and thick line) is a fit to the predictions of (optically 
thick) photoionisation models computed by Do pita et al.| (jl992 
Eq. 4. 1 therein). The thick portion corresponds to the metallicity 
range covered by the Dopita et al. models. 
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Fig. 28. Observed l ine ratios L(5007)/L(H a + [N n]) for PNe of 
the bulge of M 31 (|Ciardullo et al.||20"02| and the disk of M 33 
(Mag rini et al.|2000| l, as a function of M(5007), compared with 
the predictions of various hydrodynamical sequences with dif- 
ferent metallicities. Evolution proceeds always from faint mag- 
nitudes and low line ratios up to maximum 5007 A emission and 
then back again, but at higher line ratios. Foreground extinction 
is considered for the data, but not internal extinction. The arrow 
indicates the direction of shifts to be applied to the tracks if ad- 
ditional internal extinction occurs. 



10.0 



1 



o 
o 

ID 



i I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | 



I I I I I I I I I 



0.696 M s , Zqq 

0.696 M s , Z GD /3 

0.696 M s , Z GD /10 

■ - ■ 0.625 M S , Zqq 



0.625 M e ,Z GD /3 

0.625 M S ,Z GD /10 

— 0.605 M S , Zqq 

— - 0.595 M , Zqq 




-3 -2 -1 
M(Ha+[NII]) 

Fig. 29. The same data and theoretical sequences as presented 
in Fig. 28 but now plotted over M(Ha + [N n]). Objects and se- 



quences are only plotted down to M(5007) = 1, the brightness 
limit of the samples (cf. Fig.[28|. 
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Fig. 30. The ratios L(5007)/[L(3726) + L(3729)] are shown 
against L(4686)/L(H/3). The predictions from a selection of our 
simulations (the two 0.605 M0 sequences are from Paper IV) is 
compared with Local Group data taken from Richer ( 2006jref- 
eren ces therein and priv. co mm.). The Galactic bulge data are 
from Stasihska et al. ( 1998| l. Shown are only objects with well- 
determined line strengths and which are not fainter than two 
magnitudes below the [O in] PNLF cut-off. Evolution proceeds 
from left to right and then down because of recombination. In 
order to be consistent with the data, the tracks are truncated at 
two magnitudes below the PNLF cut-off. The kinks seen in two 
sequences are due to the optically thick/thin transition. 
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Fig. 31. Evolution of PN models in the V H whm-^(5007) plane. 
Shown are nebular sequences around central star models of 
0.696 and 0.625 Mq, each of which started with the same ini- 
tial model but with three different metal contents, as indicated 
in the figure legend. These sequences are supplemented by ones 
with less massive central stars of 0.605 and 0.595 Mq, both with 
hydrodynamically computed initial envelopes (cf. Table [2]). The 
abscissa gives either the [O m] magnitude M(5007) as defined by 
Jacoby|(fl989} (top) or M(5007) - M*, with M* = -4.48 (bot- 
tom). The evolution proceeds from low velocities and faint mag- 
nitudes over medium velocities and bright magnitudes towards 
high velocities and again faint magnitudes. The models are com- 
pared with a large data set of bright PNe belonging to a variety 
of galaxies (or parts of them) as indicated in the legend and ex- 
tracted from |Richerl ( |2006| Fig. 4 therein). 
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Fig. 34. V H whm(5007) vs. M(5007) for 0.696 M sequences 
(a = 2) with the AGB mass-loss rates (top) and AGB wind ve- 
locities (bottom) varied, as indicated in the legends. In the top 
panel V ag b = 15 kms -1 holds in all cases, whereas we have al- 
ways M ag b = 1 x 10~ 4 M0yr~' in the bottom panel. Our stan- 
dard 0.696 M0 sequence with M agb = 1 x 10~ 4 M0yr"'and 
Vagb - 15 kms~' serves in both panels as a reference (solid 
lines). Evolution proceeds always from low to high velocities. 
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Fig. 35. Same as in Fig. 32 but for the sequences with 0.605 
Mq (top) and 0.595 Mq (bottom) already presented in Fig. 31 
The post-AGB age ranges considered are 12 000 years for both 
sequences. 



50 



E 



t — i — i — i — i — i — i — i — r 



0.625 M e Z GD - [0111] 

0.625 M S Z GD - Hbeta 

0.625 M e Z GD /3-[OIII] 

0.625 M e Z GD /3 - Hbeta 

0.625 IVUZnn/10-IOIII] 



0.625 M 



Hbeta 




D. Schonberner et al.: The evolution of planetary nebulae. VII 

40 



37 



1 — I — I — I — I — I — I — I — I— 



M(Hbeta) 

Fig. 36. Run of V H whm(5007) and V H whm(H/3) vs. M(H/3) for the 
different 0.625 M0 sequences. Evolution proceeds in all cases 
from faint magnitudes and low velocities towards faint magni- 
tudes and higher velocities. 
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Fig. 37. V H whm(5007) vs. V H whm(H/3) for four 0.625 M and 
two 0.605 sequences with different metallicities as indi- 
cated in the legend. Evolution proceeds always from low towards 
higher velocities. Plotted are only the evolutionary phases with 
T e « < r max . The models are compared with recent HWHM ve- 
locity measurements of Milky Way bulge objects, either with 
apparently unsplit line profiles (filled symbols) or with split 
profiles (open symbols), provided to us by M. Richer (see 



also Richer et al. 2008 


2009 


2010a), and assuming that the 


halfwidth velocities of H/3 anc 


Ha are the same. Only objects 



with 5007 A line strengths greater than Hfi are considered. 
Accordingly, thin line sections at low velocities indicate models 
for which L(5007) < L(Hfi). The model profiles are broadened 
by a Gaussian with FWHM of 11 kms -1 in order to comply 
with the spectral resolution of the observations. Crossed sym- 
bols mark objects which are known to posses [WC]-type central 
stars (see Gorny et al.|2004| l. 
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Fig. 38. Similar as in Fig. 
Vhwhm(H/0 vs. M(H/3) predicted from model sequences with 
different central stars and the same metallicity. Evolution pro- 
ceeds again from faint magnitudes and low velocities towards 
faint magnitudes and higher velocities. These theoretical pre- 
dictions are again compared with the Galactic bulge sample 
of |Richer efaL] p009l |2010a| l already used in Fig. [37] Plotted 
are only objects with unsplit lines, and thus the region con- 
taining higher velocities at low brightness is devoid of objects. 
Each object (except for a few cases) is plotted twice, with its 
Vhwhm(5007) (filled symbols) and Vhwhm(H/3) (= Vhwhm (Hot)) 
(open symbols). We assumed a distance of 7.5 kpc for all ob- 
jects. The theoretical profiles are broadened by a Gaussian with 
FWHM of 1 1 km s -1 in order to resemble the observations. 
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Fig. 40. Evolution of R out / Vhwhm (top) and 
Vhwio%m/0.911 VWhm (bottom) vs. absolute 5007 A brightness, 
M(5007), for different nebular model sequences. In both panels, 
evolution proceeds from faint magnitudes at the upper right 
towards faint magnitudes at the lower right. 



